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I.  INTRODUCTION 


Microwave  and  high-speed  very  large  .scale  integrated  (VLSI)  circuits  are  lim¬ 
ited  by  the  propagation  characteristics  of  the  on-chip  interconnections  [Ref.  1].  The 
effects  of  signal  dispersion  and  loss  become  even  more  predominant  as  devices  and 
circuits  are  scaled  to  smaller  dimensions.  At  microwave  frequencies,  the  intercon¬ 
nections  between  elements  on  a  dielectric  substrate,  such  as  silicon  (Si)  or  gallium 
arsenide  (GaAs)  where  considerable  wavelength  reductions  occur,  must  be  treated  as 
planar  waveguide  structures.  The  analysis  and  design  of  circuits  consisting  of  these 
guided  wave  structures  are  facilitated  by  the  use  of  equivalent  circuits.  Accurate 
high  frequency  characterization  of  picosecond  pulse  propagation  on  these  structures 
including  dispersion  and  losses  requires  extensive  numerical  techniques.  The  em¬ 
pirical  equations  that  have  been  derived  for  the  propagation  parameters  from  these 
numerical  results  do  not  lead  to  an  equivalent  circuit  model  with  realizable  series  and 
shunt  branches.  Furthermore,  reliable  and  accurate  empirical  models  for  many  useful 
structures  such  as  microstrips  on  dielectric  substrate  used  in  monolithic  microwave 
integrated  circuits  (MMIC’s)  are  not  available  [Ref.  2:  pp.  256-2G2]. 

Knowledge  of  the  transient  signal  behavior  on  microstrip  transmission  lines  is  es¬ 
sential  for  the  design  of  MMIC’s  at  high  switching  speeds  or  high  frequencies.  High¬ 
speed  time-domain  measurements  must  be  used  to  properly  understand  and  mode! 
the  transient  response.  The  microwave  techniques  available  for  picosecond  pulse  prop¬ 
agation  characterization  are  scattering  parameter  measurements,  high-speed  sam¬ 
pling  oscilloscope  measurements,  and  picosecond  photoconductor  measurements.  The 
scattering  or  s-parameter  measurement  method  is  a  small-signal  frequency-domain 
technique  which  is  widely  used  by  microwave  network  analyzers.  High-speed  sam¬ 
pling  oscilloscope  measurements  can  characterize  low-level  signals  and  have  a  higher 
signal-to-noise  (S/N)  ratio  than  the  s-parameter  method.  Picosecond  photoconduc¬ 
tor  measurements  are  a  new  high-frequency  measurement  technique  made  possible  by 
advances  in  laser  technology.  The  generation  of  optical  pulses  with  sub-picosecond 
durat  ion  are  shorter  than  those  that  can  be  generated  solely  by  electronic  means  [Ref. 
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3:  pp.  117-124].  An  optoelectronic  transducer  such  as  a  photoconductor  permits  the 
conversion  of  these  sub-picosecond  optical  pulses  to  picosecond  electrical  signals.  The 
significant  advantage  of  this  technique  is  the  photoconductors  can  be  integrated  on 
the  substrate  material  to  facilitate  very  high-speed  measurements  with  extremely  high 
S/N  and  sensitivity  to  microvolt  signal  levels.  The  disadvantage  of  the  s-parameter 
and  sampling  oscilloscope  measurement  methods  is  they  suffer  from  poor  connections 
in  terms  of  high-frequency  signal  transmission  on  the  substrate. 

There  is  an  urgent  need  for  the  development  of  computer-aided  design  (CAD) 
models  tailored  specifically  to  the  special  demands  of  MMIC  technology.  MMIC’s  are 
not  readily  tunable  after  manufacture,  and  the  time  and  costs  involved  in  a  design 
are  high.  In  addition  to  the  development  of  new  models  for  CAD  design,  a  diagnostic 
testing  procedure  that  permits  on-wafer  characterization  of  MMIC’s  before  dicing  the 
wafer  into  individual  chips  is  highly  desirable  [Ref.  4:  p.14]. 

It  will  be  useful  to  formulate  a  digital  signal  processing  framework  to  solve  this 
modeling  problem.  The  motivation  for  this  approach  is  simple.  The  availability  of  in¬ 
put  and  output  time-domain  measurements  invites  the  transition  from  an  equivalent 
circuit  model  to  a  parametric  model.  For  our  purposes,  we  will  define  a  parametric 
model  as  a  discrete  mathematical  description  of  the  actual  propagation  mechanism  of 
the  system  in  terms  of  specific  model  parameters  or  coefficients.  Three  particularly 
important  discrete-time  parametric  models  are  the  moving-average  (MA),  the  autore¬ 
gressive  (AR),  and  their  combination,  the  autoregressive  moving-average  (ARMA) 
model.  These  models  are  described  as  linear  transformations  of  respective  input- 
output  time  series. 

The  model- building  process  is  determined  from  a  priori  (structural)  knowledge 
and  a  posteriori  (measurement)  knowledge  of  the  system.  The  basic  approach  taken 
in  this  dissertation,  with  respect  to  the  modeling  of  microwave  integrated  circuit 
(MIC)  interconnections,  is  to  represent  them  as  linear  time-invariant  systems.  This 
allows  linear  algorithms  to  be  applied  to  the  solution  of  the  modeling  and  parameter 
estimation  problem.  The  structural  knowledge  of  the  microstrip  line  permits  the  use 
of  empirical  and  analytical  expressions  to  formulate  lumped-distributed  equivalent 
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circuit  models  which  can  characterize  the  system.  Rational  network  functions  are 
evaluated  from  the  circuit  models  and  discrete  linear  transformations  are  employed 
to  produce  the  parametric  models. 

Another  main  emphasis  of  this  research  is  on  the  system  identification  aspect. 
Identification  is  defined  by  Zadeh  [Ref.  5:  pp.  856-865]  as:  "the  detcrrninatwu  on  the 
basis  of  input  and  output,  of  a  system  within  a  specified  class  of  systems,  to  which  the 
system  under  test  is  equivalent.’'  Consequently,  the  identification  problem  is  reduced 
to  that  of  model  parameter  estimation  based  upon  measurement  knowledge. 

A.  HISTORICAL  BACKGROUND 

Numerous  authors  have  focused  their  research  on  various  computational 
approaches  to  analyze  and  model  microstrip  on  different  substrates  [Refs.  6,7].  Ac¬ 
curate  and  complete  analysis,  however,  requires  elaborate  mathematical  models  and 
time  intensive  numerical  techniques.  Quasi-TEM  mode  calculations  combined  with 
frequency-dependent  expressions  can  achieve  accuracies  within  one  percent  of  the  full- 
wave  analysis  [Ref.  8].  Computer-aided  design  (CAD)  programs  are  available  with 
the  capability  of  both  synthesis  and  analysis  of  microstrip  lines.  The  extensive  de¬ 
velopment  of  closed-form  expressions  for  the  frequency-dependent  effective  microstrip 
permittivity  «c//(/),  has  produced  rapid  computation  of  the  microstrip  parameters. 
However,  the  empirical  equations  that  have  been  derived  for  the  frequency-dependent 
propagation  terms  do  not  lead  to  an  equivalent  circuit  model  with  realizable  series  and 
shunt  circuit  impedances.  It  is  these  equivalent  impedances  that  relate  to  the  physi¬ 
cal  microstrip  properties.  Therefore,  the  effective  microstrip  permittivity  of  different 
dielectric  substrates  is  critical  for  the  synthesis  and  validation  of  IC  interconnections 
models. 

Recently,  optical  techniques  have  been  successfully  used  in  the  characterization 
of  microwave  devices  and  integrated  circuits.  Frequency-domain  measurements  have 
been  performed  using  electro-optic  probing  of  a  microstrip  line.  In  this  work,  the 
microwave  signal  is  launched  onto  the  circuit  using  coplanar  waveguide  (CPW)  con¬ 
tacting  probes.  Hung  et  al.  [Ref.  9]  have  developed  an  on-wafer  GaAs  MMIC 


measurement  system  using  a  picosecond  pulse  sampling  technique  proposed  by  Aus- 
ton  [Ref.  10].  This  optoelectronic  characterization  has  been  demonstrated  to  achieve 
a  broad-band  frequency  response  for  both  the  magnitude  and  the  phase  of  a  K„-band 
(27-40  GHz)  MMIC. 

The  problem  of  system  identification  is  the  determination  of  a  mathematical 
model  that  provides  an  optimum  characterization  for  a  system  or  process.  Usually 
our  knowledge  of  the  system  is  limited  by  observable  input-output  measurements. 
The  inverse  scattering  problem  as  applied  to  electromagnetic  systems  has  close  con¬ 
nections  to  several  signal  processing  concepts  such  as  the  design  of  digital  filters,  the 
development  of  linear  prediction  algorithms  and  their  lattice  filter  implementations. 
In  particular,  transmission-line  systems  on  planar  dielectric  materials  can  be  inter¬ 
preted  as  layered,  wave  scattering  structures.  Two  successful  analogous  applications 
have  been  in  speech  and  seismic  signal  processing.  In  speech,  the  vocal  tract  has  been 
modeled  as  an  acoustic  tube  of  varying  cross-sectional  area  [Ref.  11].  The  inverse, 
or  model  identification,  problem  is  to  determine  the  medium  properties  from  its  re¬ 
flection  response  measured  at  an  observable  boundary,  to  some  incident  input  signal. 
This  also  relates  to  the  seismic  problem.  Here,  an  impulsive  input  to  the  earth,  such 
as  an  explosion,  will  generate  seismic  waves  propagating  downwards.  Reflections  are 
produced  as  the  wave  encounters  respective  earth  layers.  Prom  the  resulting  reflected 
waves,  the  layered  structure  beneath  the  surface  can  be  identified. 

Bruckstein  and  Kailath  [Ref.  12]  have  specified  conditions  on  the  layered  structure 
that  permits  the  use  of  recursive  parameter  estimation  algorithms.  Their  first  algo¬ 
rithm  takes  scattering  data  and  processes  it  to  identify  a  unique  layer  and  then,  at  each 
iteration,  replaces  the  data  by  a  set  of  “synthetic”  scattering  data.  This  procedure 
is  called  the  layer-peeling  method.  The  second  approach,  called  the  layer-adjoining 
method,  propagates  the  original  scattered  data  through  the  previously  identified  layer 
to  determine  information  about  future  layers.  The  layer-adjoining  method  is  the 
Levinson-Durbin  alogrithm  for  solving  the  Toeplitz  normal  equations.  For  lossless 
structures,  the  layer-peeling  method  is  the  Schur  algorithm.  The  primary  difference 
between  the  layer-adjoining  and  the  layer-peeling  process  is  that  the  latter  method 
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has  avoided  the  requirement  to  compute  inner  products.  These  algorithms  will  pro¬ 
vide  background  for  the  derivation  of  model  parameter  estimation  algorithms  that 
will  be  presented  in  Chapter  V. 

Bruckstein,  Levy  and  Kailath  [Ref.  13]  have  discussed  several  classes  of  physical 
models  that  are  equivalent  descriptions  of  a  lossless  scattering  media.  They  presented 
how  a  lumped  circuit  model  for  a  uniform  lossless  transmission  line,  described  by  par¬ 
tial  differential  equations,  can  be  formulated  into  a  discretized  wave  scattering  layered 
model.  This  model  illustrates  how  propagating  waves  through  a  lossless  transmission 
line  can  be  described  as  a  multi-layered  medium.  Each  layer  is  characterized  by  its 
impedance  function  which  directly  relates  to  a  layer  reflection  coefficient.  In  Chap¬ 
ter  IV,  multiple  microstrip  discontinuities  are  compared  to  the  multi-layered  model; 
however,  loss  and  dispersion  are  considered. 

B.  OBJECTIVES  OF  THE  RESEARCH 

The  first  objective  of  this  research  is  to  develop  an  equivalent  circuit  model  for 
a  dispersive  lossy  microstrip  transmission  line  which  is  compatible  with  the  standard 
circuit  analysis  and  design  techniques,  including  computer-aided  design  (CAD)  tools. 
The  proposed  circuit  model  will  then  be  extended  to  describe  abrupt  width  dimen¬ 
sion  changes,  called  impedance  discontinuities,  of  the  microstrip  line.  The  transient 
analysis  of  these  equivalent  circuits  will  simulate  the  picosecond  pulse  propagation 
on  dispersive  lossy  microstrip  lines  and  the  effects  of  impedance  discontinuities  can 
be  modeled.  The  circuit  models  will  be  verified  using  photoconductor  measurements 
from  IC  interconnections  and  experimental  results  from  fabricated  microstrip  test 
structures. 

The  second  objective  of  the  thesis  is  to  derive  ARMA  and  AR  parametric  mod¬ 
els  for  three  microstrip  test  structures.  These  include  1)  an  impedance  matched 
microstrip  transmission  line,  2)  a  cascaded  microstrip  step  discontinuity,  and  3)  a 
multi-section  microstrip  step  discontinuity.  The  impedance  matched  microstrip  line 
describes  the  typical  MIC  interconnection  [Ref.  3:  pp.  18-19].  As  discussed  eariler, 
the  formulation  of  these  parametric  models  provides  an  opport  unity  to  exploit  existing 
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parameter  estimation  algorithms,  as  well  as  develop  new  algorithms  to  characterize 
transient  signal  propagation  on  MIC’s. 

The  last  objective  is  the  development  of  parameter  estimation  algorithms.  This 
research  work  will  focus  on  the  estimation  of  both  ARMA  and  AR  model  parameters 
from  die  impulse  response  of  an  electrically  short  microstrip  line  section.  Addition¬ 
ally,  the  system  identification  of  the  multi-section  microstrip  structure  will  also  be 
investigated.  The  detection  of  each  section  impulse  response  will  be  complicated 
by  pulse  dispersion,  propagation  loss,  and  multiple  reflections  due  to  the  impedance 
discontinuities.  A  new  layer-probing  algorithm  will  be  presented  to  overcome  these 
difficulties.  Computer  simulations  of  the  algorithms  are  performed  to  validate  their 
performance. 

C.  ORGANIZATION  OF  THE  REPORT 

Chapter  II  derives  a  lumped-distributed  equivalent  circuit  that  includes  the  ef¬ 
fects  of  microstrip  loss  and  dispersion.  The  bandwidth  of  the  propagating  pulse  will 
establish  a  maximum  wavelength  of  interest.  The  physical  length  of  the  IC  intercon¬ 
nection  being  considered  will  be  very  much  less  than  a  quarter- wavelength.  Under 
this  condition  the  lumped-element  circuit  approximation  is  used  to  characterize  the 
microstrip  line  section. 

Chapter  III  presents  the  parametric  models  for  the  impedance  matched  quarter- 
wavelength  microstrip  transmission  line.  An  ARMA  parametric  model  is  derived 
from  the  network  function  of  the  equivalent  circuit  model.  The  ARMA  model  coef¬ 
ficients  are  shown  to  be  directly  related  to  the  equivalent  lumped  capacitance  and 
inductance  of  the  circuit  model.  An  AR  parametric  model  will  also  be  presented. 
Relationships  are  given  to  evaluate  the  effective  microstrip  permittivity  directly  from 
the  ARMA/AR  model  parameters. 

In  Chapter  IV,  microstrip  discontinuities  will  be  considered.  Equivalent  circuit 
models  are  derived  for  a  cascaded  microstrip  step  discontinuity.  Similarly,  ARMA  and 
AR  parametric  models  are  developed  using  discrete  transformations  of  the  equivalent 
circuit  network  functions. 
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Parameter  estimation  algorithms  are  the  emphasis  in  Chapter  V.  Several  esti¬ 
mation  algorithms  are  developed  for  both  deterministic  and  stochastic  data.  Three 
deterministic  based  algorithms  are  presented.  A  weighted  least  squares  (WLS)  algo¬ 
rithm  is  derived  to  solve  the  ARMA  model  parameters  from  a  finite-length  impulse 
response.  An  alternative  technique  will  analytically  approximate  a  transfer  function 
of  the  microstrip  section  using  measured  input-output  rise  and  delay  times  of  a  tran¬ 
sient  pulse.  Finally,  AR  model  parameters  are  estimated  from  the  impulse  response 
using  the  Schur  algorithm  [Ref.  14].  The  stochastic  algorithms  will  assume  that  a 
white  noise  source  is  applied  to  the  input  of  the  microstrip  structure.  The  Schur 
algorithm  will  estimate  the  AR  model  parameters  from  the  output  data.  When  input 
and  output  random  data  are  available,  an  ARMA  parameter  estimation  algorithm 
based  on  a  generalized  Mullis- Roberts  (M-R)  criterion  [Ref.  15]  is  employed.  How¬ 
ever,  if  only  the  output  data  can  be  measured,  a  modified  two-stage  least  squares 
algorithm  is  presented  to  estimate  a  second-order  ARMA  model.  Finally,  new  sys¬ 
tem  identification  algorithms  for  the  multi-section  microstrip  step  discontinuity  are 
presented. 

Chapter  VI  presents  simulation  and  experimental  results  Waveform  compar¬ 
isons  are  made  between  measured  data,  equivalent  circuit  simulations,  and  paramet¬ 
ric  mode!  simulations  for  each  microstrip  test  structure.  The  performance  of  the 
parameter  estimation  algorithms  will  be  investigated. 

Chapter  VII  is  a  summary  of  the  significant  contributions  presented  in  this  dis¬ 
sertation.  It  draws  conclusions  from  the  results  and  proposes  some  important  future 
directions  for  this  research. 

Three  appendices  are  included.  Appendix  A  contains  an  alternate  proof  of  the 
Schur  algorithm.  Appendix  B  contains  the  equivalent  circuit  model  listings  that 
are  used  by  the  transient  analysis  program.  Appendix  C  contains  listings  of  the 
FORTRAN  programs  used  in  the  simulations  presented  in  this  report. 
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II.  MICROSTRIP  CHARACTERIZATION 


In  this  chapter  equivalent  circuit  models  are  derived  for  a  lossy  dispersive  mi¬ 
crostrip  transmission  line.  The  propagation  characteristics  of  the  transmission  line 
are  modeled  by  lumped-distributed  equivalent  circuit  models.  In  order  to  derive  the 
expressions  for  the  equivalent  lumped  circuit  elements,  the  models  are  defined  for 
a  maximum  frequency  of  interest  which  is  restricted  by  the  physical  length  of  the 
microstrip  line  and  the  dielectric  constant  of  the  substrate.  The  empirical  equations 
described  in  the  literature  for  the  propagation  parameters  do  not  usually  lead  to  an 
equivalent  circuit  model  [Ref.  16]. 

The  equivalent  circuit  models  presented  are  compatible  with  standard  circuit 
analysis  and  design  techniques  including  the  use  of  computer-aided  design  tools  such 
as  PSPICE  [Ref.  17].  The  proposed  models  will  include  the  effects  of  dispersion 
and  loss  at  microwave  frequencies. 

A.  IC  MICROSTRIP  INTERCONNECTIONS 

Integrated  circuit  interconnections  can  be  described  by  microstrip  transmission 
lines  because  their  geometries  are  similar.  The  abrupt  dielectric  interface  shown 
by  the  open  microstrip  geometry  in  Figure  2.1  makes  it  incapable  of  supporting 
a  single  model  of  propagation.  However,  microstrip  propagates  the  bulk  of  its 
energy  in  a  field  distribution  which  approximates  the  transverse  electromagnetic 
(TEM)  mode  and  is  usually  referred  to  as  the  quasi-TEM  or  quasi-static  mode 
[Ref.  18].  Several  computational  approaches  are  available  utilizing  quasi-TEM  mode 
calculations  combined  with  closed-form  frequency-dependent  expressions  [Ref.  6]. 
It  is  this  latter  approach  which  is  considered  in  this  work. 
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The  characteristic  impedance  of  a  TEM  transmission  line  is  described  by 


Z0  = 


(2.1) 


where  L  and  C  are  the  equivalent  lumped  inductance  and  capacitance  of  the  mi¬ 
crostrip  section,  respectively.  With  microstrip  geometries,  the  same  type  of  di¬ 
electric  substrate  is  used  below  the  conductor.  However,  there  is  air  above  the 
conductor. 


Ground  Plane  (Conducting) 


Figure  2.1  Generalized  open  microstrip  geometry 

The  dielectric  constant  used  in  the  design  must  take  into  consideration  the  dielectric 
constant  of  air  (er  =  1)  and  that  of  the  substrate  material.  The  effective  microstrip 
permittivity  teff  will  be  .ntroduced.  This  quantity  is  unique  to  mixed-dielectric 
transmission  line  systems  and  it  provides  a  useful  link  between  different  wavelengths, 
impedances,  and  propagating  velocities. 
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An  expression  for  the  static-TEM  effective  microstrip  permittivity  has  been  calcu¬ 
lated  by  Owens  [Ref.  19]: 


(2.2) 


where  er  is  the  relative  dielectric  constant  of  the  substrate.  Using  the  effective  mi- 
crostrip  permittivity,  the  characteristic  impedance  at  TEM  frequencies  is  calculated 
by 
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where  w  and  h  are  microstrip  conductor  width  and  substrate  thickness,  respectively 
[Ref.  20].  A  correction  factor  is  applied  to  account  for  the  fringing  fields  associated 
with  a  finite  conductor  thickness.  Bahl  et  al.  [Ref.  21]  have  introduced  an  effective 
width  parameter  into  (2.3)  in  order  to  improve  the  Z0  calculation.  The  effective 
width  we ,  replaces  the  w  in  equation  (2.3) 
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where  t  is  the  conductor  thickness. 


B.  FREQUENCY-DEPENDENT  MICROSTRIP  MODELS 

The  quasi-TEM  analysis  of  microstrip  transmission  lines  mentioned  above  be¬ 
gins  to  lose  accuracy  at  high  microwave  frequencies.  The  characteristic  impedance 
and  the  effective  microstrip  permittivity  are  dependent  upon  the  dielectric  thickness- 
to-guide  wavelength  ratio,  2nh/\g  [Ref.  22].  In  addition,  different  loss  mechanisms 
become  important  and  the  attenuation  function  is  also  frequency  dependent.  This 
frequency-dependency  of  the  microstrip  is  caused  by  a  hybrid  mode  of  propagation 
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that  describes  a  coupled  version  of  both  transverse  electric  (TE)  and  transverse 
magnetic  (TM)  modes. 

The  propagation  group  velocity  of  a  signal  also  depends  on  the  frequency- 
dependent  effective  microstrip  permittivity  as 


Mf)  = 


c 

V^IKT) 


(2.5) 


where  c  is  the  velocity  of  light  in  a  vacuum.  Pulses  which  have  a  spectral  component 
above  the  TEM  mode  frequency  regime  will  be  dispersed  because  the  higher  har¬ 
monics  of  the  pulse  will  travel  at  a  slower  phase  velocity  than  the  lower  harmonics 
[Ref.  IS].  Therefore,  the  phase  constant  ( 3 )  is  a  non-linear  function  of  frequency, 
which  leads  to  the  phenomena  of  dispersion. 

Several  methods  for  evaluating  the  frequency-dependent  effective  microstrip  per¬ 
mittivity  are  available  [Refs.  23.  24].  However.  Yamashita  et  al.  [Ref.  7]  have 
derived  an  expression  by  curve  fitting  the  data  obtained  from  a  full-wave  analysis 
as 
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If  the  substrate  material  is  lossless,  the  phase  constant  becomes 
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However,  finite  resistivity  of  the  conductor  and  finite  conductivity  of  the  substrate 
introduce  attenuation.  At  low  microwave  frequencies  the  conductor  loss  factor  oc 
is  given  by  [Ref.lS:  p.  90] 
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where  /  is  in  gigahertz.  For  higher  frequencies,  Pucei  et  al.  [Ref.  25]  have  calculated 
the  conductor  loss  as 
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and  Rs  is  the  surface  skin  resistance  in  ohms  given  by 

Rs  = 
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where  //  is  the  conductor  permeability  and  crc  is  the  conductor  conductivity.  An 
expression  for  dielectric  loss  qj  has  been  derived  by  Hammerstad  and  Bekkadal 
[Ref.  26]  as 
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where  tan  6  is  the  less  tangent  of  the  substrate.  For  lossy  semiconductor  substrates, 
the  general  attenuation  and  phase  constants  for  the  microstrip  are  derived  from  an 
analysis  of  a  transmission  line  section  modeled  by  series  and  shunt  resistances  to 
represent  the  conductor  and  dielectric  loss  mechanisms  [Ref.  27]  as 
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C.  PICOSECOND  PULSE  PROPAGATION 

One  of  the  research  objectives  is  to  accurately  characterize  picosecond  pulses 
propagating  on  microwave  IC  interconnections.  Researchers  have  studied  the  ef¬ 
fects  of  dispersion  on  picosecond  pulses  by  performing  computer  simulations  in  the 
frequency  domain  [Refs.  28.29].  An  input  time-domain  signal  v(t,  0)  is  transformed 
using  Fast  Fourier  transform  (FFT)  techniques  to  produce  the  input  frequency  spec¬ 
trum.  F  jr(LO)  j.  A  frequency- dependent  propagation  function  7 (/)  acts  upon  each 
spectral  component  as 

F{r(L0)}t-7^  1 

where  the  propagation  function  is  q(/)  =  o(/)  +  j3(f)  and  /  is  the  propagation 
distance  along  the  microstrip.  The  distorted  temporal  pulse  is  determined  using  an 
inverse  Fourier  transform  as 

v(t.l)  =  F_1{F{t>(<,0)}c“lQW  iy 

The  approach  tahen  in  this  work  is  similar,  but  the  exponential  propagation  term 
is  replaced  by  a  discrete  parametric  model  of  the  microstrip  section.  Specifically, 
the  phase  term  (3(f)  is  modeled  by  a  digital  transfer  function  H(z),  while  a  discrete 
loss  multiplier  describes  the  attenuation  constant  <*(/)  of  the  microstrip  line  section 
at  a  maximum  frequency  of  interest.  Chapter  III  will  introduce  the  parametric 
models  for  several  microstrip  transmission  line  structures. 
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D.  DERIVATION  OF  EQUIVALENT  CIRCUIT  MODELS 

A  microstrip  transmission  line  is  approximated  by  cascading  several  electrically 
short  sections  [Ref.  30].  A  single  section  of  microstrip  line  having  a  characteristic 
impedance  Zq  and  terminated  in  a  matched  load  impedance  is  shown  in  Figure 
2.2a.  The  load  impedance  Z(n+i),  is  the  equivalent  input  impedance  of  the  (n  -f  1)- 
th  section.  The  normalized  input  impedance  to  the  n-th  section  is  derived  [Ref.  30] 
as 


Z„  _  Zn+ 1  cosh(7/)  +  Zq  sinh(7/) 
Zq  Zn+ 1  sinh(7/)  +  Zq  cosh(7/) 


where  the  propagation  function,  7  =  0  +  j/?,  describes  the  attenuation  a,  and  phase 
shift  /?  of  the  n-th  section.  Equation  (2.12)  is  rewritten  as 


7  +  tanh(7’/) 

__  Zj q _ 

1  4-  tanh(7 /) 

Zq 

The  hyperbolic  function,  tanh(7/),  is  expanded  to 


(2.13) 


tanh(7/)  = 


tanh(a/)  -f  j  tan (/?/) 
1  +  j  tanh(aZ)  tan (/?/) 


Substituting  (2.14)  into  (2.13)  yields 

z.  (^r) +  tanh(ai) 

1  +  ~  ^  tanh(o/) 

where 

/Zn+iV  = 

^  Z“  ‘  +  >  ““(W 


(2.14) 


(2.15a) 


(2.156) 


A  physical  interpretation  of  (2.15)  is  illustrated  in  Figure  2.2b.  The  original  n- 
th  section  is  equivalent  to  three  sub-sections  of  the  same  characteristic  impedance 
cascaded  together.  The  first  and  third  sub-sections  provides  attenuation  but  no 
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phase  shift,  and  are  modeled  by  purely  resistive  attenuator  pads.  The  total  loss  of 
the  microstrip  line  is  divided  equally  between  each  7r-attenuator.  The  center  sub¬ 
section  has  a  pure  imaginary  propagation  function  and  no  attenuation.  This  type 
of  reactive  network  corresponds  to  a  lossless  transmission  line  having  a  phase  shift 
equal  to  the  original  section. 

In  classical  filter  theory  the  reactive  lossless,  two-port  network  has  been  de¬ 
scribed  as  a  constant-k  lowpass  filter.  Two  variations  of  this  filter  network  will  now 
be  presented.  A  reactive  7r-network  is  used  to  derive  the  lumped-parameter 


Z 

o 


(a) 


0>) 

Figure  2.2  Microstrip  line  section  equivalent  circuit 
equivalent  circuit  model  and  a  reactive  T -network  is  used  to  derive  a  lumped- 
distributed  equivalent  circuit  model. 
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1.  Lumped-Parameter  Equivalent  Reactive  II-Network 

The  lossless  reactive  network  of  Figure  2.2b  is  modeled  by  an  equivalent 
7r-network.  The  complex  impedances  are 


Zi  = 


Z2  =  jZo  sin(/?/). 


(2.16a) 


(2.166) 


Equations  (2.16a)  and  (2.16b)  represent  a  capacitive  and  an  inductive  lumped  re¬ 
active  impedance,  respectively.  Solving  for  the  susceptive  {Be)  and  reactive  (Xi) 
lumped  elements  yields 


Be  =  J^tan  (y)  (2.17a) 

Xl  =  Zq  sin{/3l).  (2.176) 


Assuming  a  very  short  section  length  (/  <<  A5/4),  the  equivalent  shunt  capacitance 
and  series  inductance  can  be  approximated  as  follows  [Ref.  IS] 


Cr 


l 

2Z0  /  \g 


(2.  lSa ) 


L* 


Z0l 

\f 


(2.1S6) 


The  guide  wavelength  \g  is  expressed  in  terms  of  the  frequency- dependent  microstrip 
permittivity  ee//(/)  as 


JyfciZf) 


(2.19) 


where  /  is  the  highest  frequency  of  interest  of  the  propagating  signal  [Ref.  18: 
p.  70].  The  resultant  lumped-parameter  equivalent  circuit  for  the  microstrip  line 
section  is  shown  in  Figure  2.3. 
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2.  Lumped-Distributed  Equivalent  Reactive  T-Network 

The  reactive  7r-network  of  Figure  2.3  can  be  replaced  by  an  equivalent  T- 
network  as  shown  in  Figure  2.4a.  The  series  inductive  and  the  shunt  capacitive 
impedances  are  then  given  by 


Z\  =  jZa  tan  (  y 


Z2  = 


Z0 


j  sin(pl) ' 

Solving  for  the  reactive  and  susceptive  lumped  elements  of  Figure  2.4  yields 

ttI 


Xi  =  2Zo  tan  (  — 

Ar 


BC  =  ~  sin 
Z  o 


9 

2tt  l 
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(2.20a) 

(2.206) 

(2.21a) 

(2.216) 


Figure  2.3  Lumped- parameter  equivalent  circuit  model 


Assuming  a  section  length  that  is  less  than  a  quarter  wavelength,  the  lumped  ele¬ 
ments  can  be  approximated  as  follows  [Ref.  18]: 


Lt 
Ct  ' 


Z0l 

f*9 

l 

ZofK 


(2.22  a) 
(2.226) 


Next  an  equivalency  between  the  series  lumped  inductance  (Lj)  and  a  lossless 
transmission!  line  section  is  derived.  Referring  to  Figure  2.4b.  we  define  equivalent 
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impedances,  Zeq ,  with  respect  to  the  shunt  capacitance  Cj\  The  load  impedance 
Zi  is  equivalent  to  the  characteristic  impedance  Zq.  Therefore,  the  equivalent 
impedance  yields 

Zeq  =  Zi  +  juiLj-  (2.23a) 


Using  (2.15),  the  n-th  section  input  impedance  for  a  lossless  (a  =  0),  high  impedance 
transmission-line  section  is 


z  _  Zn+i  +jZ0  tan (fil) 
z0 

Assmning  |/5/|  <C  1,  the  input  impedance  can  be  approximated  by 


(2.236) 


Z„  ~  Zn+ 1  +jZa((3l). 


(2.23c) 


Equation  (2.23c)  is  now  equated  to  the  equivalent  input  impedance  Zeq. 


Figure  2.4(a)  Equivalent  T-network  equivalent  circuit 

Assuming  a  small  section  length  and  taking  the  imaginary  terms  in  (2.23a)  and 
(2.23c)  into  account  yields 

Lt  *  Z0  (~\  =  Z0  U  (2.24) 

where  tj  is  the  propagation  delay  through  the  lossless  transmission-line  section. 


IS 


Each  lumped  inductance  Lj  can  now  be  replaced  by  an  equivalent  distributed  trans¬ 
mission  line  as  shown  in  Figure  2.5. 


Figure  2.4(b)  Equivalent  T-network  equivalent  circuit  (continued) 


Figure  2.5  Lumped-distributed  equivalent  circuit  model 
The  accuracy  of  the  circuit  model  is  dependent  on  the  microstrip  line  length 
/,  and  the  maximum  frequency  bandwidth  u ;mox  of  the  propagating  transient  sig¬ 
nal.  The  equivalent  lumped-distributed  circuit  presented  in  this  section  is  valid  for 
electrically  short,  high  impedance  transmission  lines  which  satisfy  the  conditions 
id  <  1  /umM  and  \(3l\  <  1. 
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E.  ESTIMATION  OF  CHARACTERISTIC  IMPEDANCE 

Picosecond  pulses  propagating  on  the  IC  interconnections  are  described  as  finite 
energy  signals  and  their  maximum  cumulative  energy  is  defined  by  [Ref.  31:  p.  35] 

+oo 

Em„,  =  £  |1>+M|2  (2.25) 

n=-oo 


where  v+[l,  n]  is  the  sampled  right-propagating  voltage  measured  at  the  /- th  length 
along  the  microstrip.  To  discuss  pulse  propagation  on  a  lossy  microstrip  line,  we  first 
obtain  the  transmission- line  equations  by  applying  Kirchoff’s  voltage  and  current 
laws  to  the  alternative  distributed  equivalent  circuit  as  shown  in  Figure  2.6.  When 
low-loss,  high  frequency  conditions  are  assumed,  the  general  solution  for  the  voltage 
and  current  at  the  /-th  length  is  solved  as  [Ref.  32:  pp.  437-445] 


v(l,t) 


v+(‘-^)+v'(l+i) 


(2.26a) 

(2.266) 


where  the  superscripts  (+)  and  (-)  denote  the  right-  and  left-propagating  voltage 
waves,  respectively.  When  we  compare  these  equations,  the  right-propagating  cur¬ 
rent  and  voltage  are  given  by 


(2.27  a) 
(2.276) 


where  vg  is  the  group  velocity.  The  instantaneous  power  associated  with  the  right- 
propagating  wave  is 


p+(l,t)  -  v+(l,t)  i+{l,t) 


U( 

L  \ 

V»J  J 

2o 


(2.28a) 

(2.286) 


20 


Since  the  pulse  rise  time  is  very  fast  (<  10  ps),  the  derivative  of  the  cumulative 
energy  curve  (i.e.  the  instantaneous  power)  is  approximately  constant.  Therefore, 
the  maximum  power  of  a  causal,  transient  pulse  is  given  by 

1  N 

-Pm.,  «  (2.29) 

n=0 

where  N  is  the  number  of  voltage  sample  values.  By  equating  (2.28b)  and  (2.29). 
we  have  the  characteristic  impedance  given  by 

Z0  =  l-^n]?ea*r  (2.30) 

max 

where  t>+[/,  n]peajt  is  the  peak  amplitude  voltage  sample  of  the  propagating  pulse. 


i(f.t) 

o—  » - — 


v(M) 


MAC) 

-w- 


i  ((  +  A(.t) 

►  o 


£(  A/) ; 


C(  Ar)=T=  v(f+A  M) 


(AO 


(Lumped-elements  expressed  as  per  unit  length) 


Figure  2.6  Lossy  transmission  line  equivalent  circuit  (After  Ref.  32) 


In  summary,  distributed-lumped  equivalent  circuit  models  were  derived  for  an 
electrically  short  length  of  lossy  microstrip  transmission  line.  The  lumped  circuit 
elements  are  evaluated  at  the  maximum  frequency  of  a  bandlimited  input  signal. 
The  validity  of  the  equivalent  circuit  models  are  controlled  by  the  physical  length 
of  the  microstrip  line.  Therefore,  we  are  interested  in  picosecond  pulse  widths 
and  specific  IC  interconnections  lengths  which  satisfy  the  lumped-element  modeling 
criterion.  Chapter  III  will  use  the  proposed  circuit  models  as  a  foundation  to  derive 
the  parametric  models  for  an  IC  interconnection. 
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III.  PARAMETRIC  MODELING 


The  rationale  for  using  a  parametric  modeling  approach  to  approximate  a  mi¬ 
crostrip  transmission  line  is  intimately  related  to  the  identification  of  the  physical 
microstrip  properties.  This  chapter  will  present  both  autoregressive  moving- average 
(ARMA)  and  autoregressive  (AR)  models  of  the  equivalent  reactive  networks  pre¬ 
viously  derived  for  the  IC  interconnection.  A  second  objective  is  to  show  how 
the  model  parameters  are  related  to  the  effective  microstrip  permittivity  ee//\  the 
characteristic  impedance  Zq,  and  the  propagation  group  delay  rg  of  the  microstrip 
structui  e. 

A.  RATIONAL  TRANSFER  FUNCTION  MODELS 
1.  ARMA  Parametric  Model 

In  this  model,  the  discrete-time  input  sequence  x[n],  and  the  microstrip  sec¬ 
tion  output  sequence  j/[n],  are  related  by  the  linear  difference  equation  [Ref.  31] 

q  v 

y[n]  =  ^6*x[n  -  k\  -  ^aky[n  -  *].  (3.1) 

k=0  k=l 


where  bk  are  the  moving-average  (MA)  model  parameters  and  ak  are  the  autore¬ 
gressive  (AR)  model  parameters.  The  flowgraph  representation  of  (3.1)  is  shown  in 
Figure  3.1.  This  difference  equation  describes  an  ARMA  model  of  order  (p,  q)  and 
the  corresponding  system  transfer  function  is 


Harma{z)  = 


Yj  h* 

k= 0 


-k 


1  -  YakZ  k 

k= 1 


B(z) 

A(z) 


(3.2) 
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2.  AR  Parametric  Model 


When  all  the  bk  coefficients,  except  bo  =  1,  are  zero  in  the  ARM  A  model, 
equation  (3.1)  reduces  to 

y[n]  =  x[n]  -  ^  aky[n  -  fc]  (3.3) 

k=l 


and  the  corresponding  AR  or  all-pole  transfer  function  is  given  by  [Ref.  33:  p.lll] 


Har{z)  =  - ? - 

1  -  tlak2-k 


Figure  3.1  Realization  of  the  ARMA/AR  parametric  models 

In  Figure  3.1,  a  flowgraph  realization  of  this  model  is  shown  by  the  weighted  feed¬ 
back  signal  path. 
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B.  IC  INTERCONNECTION  PARAMETRIC  MODELS 


1.  Derivation  of  the  ARMA  Model 

The  network  transmission  function  for  the  reactive  7r-network  of  Figure  2.3 


yields 


H(s)  = 


Jl 


s 2  4-  +  u>q 


(3.5) 


where  u>o  =  (LXCK)~0'5  and  £  =  (2ZoCTu;o)-1 .  This  network  function  describes  a 
bandlimited  frequency  response  for  a  second-order  lowpass  filter.  Here  we  use  the 
impulse  invariant  design  method  to  transform  (3.5)  into  a  digital  transfer  function 
[Ref.  34].  Equation  (3.5)  can  be  expressed  in  the  partial-fraction  expansion  form  as 


H(s)  =  J2 


Jfc=l 


Ck 

s  -  dk 


(3.6) 


where 

Ci,  2 

and  the  complex  poles  are 


^’Q  -±*r/2 

2n/R^ 


dl,2  =  -C^o  ±  ju?o\/l  -  C2- 


Taking  the  inverse  Laplace  transform,  the  corresponding  unit  impulse  response  be¬ 
comes 

2 

h(t)  =  for  <>0  (3.7a) 

k= l 

and  its  discrete  representation  is 

2 

h\n}  =  Y.CktdknT  (3-76) 

Jfc=l 


where  T  is  the  sampling  interval. 
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It  should  be  observed  that  for  high  sampling  rates  (say,  T  =  1.25  ps)  the  digital 
filter  has  an  extremely  high  gain.  For  this  reason  (3.7b)  is  expressed  as 

2 

h[n)  =  Tj2CkedknT.  (3.8) 

i=i 


The  digital  transfer  function  is  then  obtained  by  taking  the  z-transform  of  (3.8)  as: 

Harma(‘)  =  g  (3'9) 

After  rearranging  the  terms  and  simplifying,  we  have 

Harma(z)  =  7— -ps  (31°) 

1  —  CL\Z  1  —  tt2  Z  ‘ 

where  the  filter  coefficients,  also  called  the  ARMA  model  parameters,  are: 

bl  =  pT  e-C^o  T  cos  (^-uoy/l-C2  t)  (3.11a) 

v  l  —  C2 

a\  —  2  e— cos  (3.116) 

a2  =  _e~2C^0  T.  (3.11c) 

The  filter  coefficient  6i  is  directly  linked  to  the  AR  filter  coefficients  a\  and  a2. 
Given  a2  and  the  sampling  period  T,  (3.11c)  yields  £u>g-  Substituting  this  result 
into  (3.11b),  a  value  for  u>o  \J\  —  (2  can  be  found.  Finally,  specific  values  for  ujq  and 
(  are  easily  determined  and  the  numerator  filter  coefficient  6i  can  be  calculated. 
Therefore,  the  MA  model  parameter  in  (3.10)  can  be  obtained  from  the  AR  model 
parameter  estimates.  In  Chapter  IV  we  will  take  advantage  of  this  computational 
result  in  our  discussion  of  parameter  estimation  algorithms. 
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2.  Derivation  of  the  AR  Model 


The  AR  parametric  model  is  derived  from  the  lumped-distributed  equivalent 
circuit  model  of  Figure  2.5.  A  network  reflection  function  T(s)  is  solved  as  [Ref.  35] 


=  vw-noo 

Y0(s)  +  Yl(s) 


(3.12) 


where  the  equivalent  admittances  are  given  by 

VoW  =  A 
Yds)  =  Yc(s)  +  Yds) 
=  sCt  +  ■ 


After  simpilfying  and  rearranging  (3.12),  we  have 


The  network  transmission  function  is  then  defined  as 

T(s)  =  1  +  I». 


(3.13) 


(3.14) 


By  substituting  (3.13)  into  (3.14),  a  first-order,  single  time  constant  transfer  func¬ 
tion  is  achieved: 


T(s)  = 


6 

s  -f  6 


(3.15) 


where  6  — 
response 


Taking  the  inverse  Laplace  transform  yields  the  unit  impulse 
h.AR(t)  =  6e~ ^  *  for  t  >  0.  (3.16) 


The  /-transform  of  the  corresponding  normalized  sampled  transmission  impulse 
response  is  given  by 

Har(*)  = - b - -  for  M  >  |e-r"|.  (3.17) 

1  —  e~lft  z-1 
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where  Tn  is  the  normalized  sampling  interval.  In  order  to  satisfy  the  Nvquist  sam¬ 
pling  criterion,  the  sampling  intervals  of  the  normalized  and  denormalized  trans¬ 
mission  responses  must  be  related  by  [Ref.  36:  pp.  341-342] 

Tn  =  6T .  (3.18) 


To  adjust  the  dc  gain  of  (3.17)  we  can  multiply  the  digital  transfer  function  Har{z) 
by  the  ratio  of  (3.15)  evaluated  at  s  =  0  to  (3.17)  evaluated  at  z  =  1  (  which  is  the 
dc  value).  Thus,  the  AR  transfer  function,  after  being  adjusted  for  the  dc  gain, 


becomes 


Har(z)  = 


1  -  e~ln 

—  p~Tn  r~1 


for  \z\  >  \e~ 


(3.19) 


3.  Effective  Microstrip  Permittivity 

The  effective  microstrip  permittivity  is  directly  related  to  the  ARMA  model 
parameter  02,  the  sampling  interval  T,  and  the  physical  microstrip  length  1.  Using 
(3.11c),  the  propagation  group  delay  is  solved  as  [Ref.  18:  p.  3] 


C^O  =  - 


In  02 
2  T 

(In  o2)  (/) 


(3.20a) 

(3.206) 


where  the  group  delay  of  the  reactive  network  is  rg  =  (2 LrCr)  '  ■  The  propagation 
group  delay  can  also  be  related  to  the  effective  microstrip  permittivity  by  [Ref.  18: 

p.  62] 

Tg  =  .  (3.21) 


Substituting  (3.21)  into  (3.20b)  yields  the  effective  microstrip  permittivity 


(■rff  = 


2  T  c 
(lna2)  (/) 


(3.22) 
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We  now  derive  an  expression  for  the  effective  microstrip  permittivity  from  the 
AR  model  parameter.  From  (3.19),  we  recognize  the  first-order  AR  filter  coefficient 
is 


(3.23) 


The  effective  microstrip  permittivity  is  now  solved  in  terms  of  this  AR  model  pa¬ 
rameter  and  the  propagation  group  delay: 


C  = 


Ta  = 


In  a 


2  T 


(In  a)  (/) 


(3.24a) 

(3.246) 


Substituting  (3.21)  into  (3.24b)  yields  the  effective  microstrip  permittivity 


teff 


'  2  Tc  ' 

.(In  <0(0. 


2 


(3.25) 


4.  Modeling  Propagation  Loss 

Both  the  ARMA  and  AR  digital  transfer  functions  were  derived  from  a  loss¬ 
less  equivalent  circuit.  In  order  to  properly  characterize  a  picosecond  pulse  prop¬ 
agating  on  the  IC  interconnection,  we  must  include  a  loss  mechanism  in  the  para¬ 
metric  models.  The  total  propagation  loss  (in  dB)  experienced  by  a  high-speed 
transient  pulse  will  be  defined  by  [Ref.  37:  p.  7] 


LdB  =  101og10 


,-2««  j 


(3.26) 


where  at  is  the  total  attenuation  factor  (in  nepers/m)  and  /  is  the  propagation 
distance  (in  meters)  along  the  microstrip.  The  attenuation  factor  can  be  solved  by 
combining  the  theoretical  conductor  and  dielectric  attenuation  factors  from  (2.8b) 
and  (2.10),  respectively. 
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In  summary,  this  chapter  has  presented  ARMA  and  AR  parametric  models  for  a 
lossy  dispersive  microstrip  transmission  line.  Analytical  expressions  were  derived  to 
relate  the  effective  microstrip  permittivity  to  the  model  parameters  at  a  maximum 
frequency  of  interest.  We  have  assumed  no  discontinuities  on  the  IC  interconnection 
because  of  impedance  matching.  In  Chapter  IV,  we  will  extend  the  equivalent  circuit 
and  parametric  models  for  a  cascaded  microstrip  step  discontinuity.  Chapter  V  will 
present  several  parameter  estimation  algorithms  which  will  be  used  to  approximate 
the  effective  microstrip  permittivity  of  the  IC  interconnection. 
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IV.  MICROSTRIP  DISCONTINUITY  MODELING 


Numerous  approaches  have  been  made  to  use  equivalent  circuits  to  model  the 
microstrip  discontinuities.  Whinnery  and  Jamieson  [Ref.  38]  used  Hahn’s  method 
to  match  electromagnetic-wave  solutions  across  discontinuities,  Oliner  [Ref.  39] 
used  Babinet’s  principle  to  describe  stripline  discontinuities,  and  Menzel  and  Wolff 
[Ref.  40]  calculated  the  frequency-dependent  properties  of  various  microstrip  dis¬ 
continuities. 

A  discrete  parametric  method  is  presented  in  this  chapter  to  model  the  reflec¬ 
tion  and  transmission  characteristics  of  a  cascaded  microstrip  step  discontinuity, 
assuming  the  quasi-TEM  mode  of  propagation.  Equivalent  circuits  are  developed 
using  the  shunt  capacitance  circuit  model  presented  by  Gupta  and  Gopinath  [Ref. 
41].  Both  reflection  and  transmission  network  functions  are  derived  from  the  circuit 
models.  ARMA  and  AR  digital  filters  are  described  using  discrete  transformations 
of  the  respective  network  functions.  This  method  has  the  advantage  that  complex 
microstrip  circuits  containing  discontinuities  can  now  be  modeled  by  linear  differ¬ 
ence  equations. 

A.  MICROSTRIP  DISCONTINUITY  EQUIVALENT  CIRCUITS 
1.  Single  Step  Discontinuity 

A  single  microstrip  step  discontinuity  exists  at  the  junction  of  two  microstrip 
lines  having  different  widths  and  characteristic  impedances.  This  type  of  discon¬ 
tinuity  is  seen  in  the  design  of  microwave  matching  transformers,  couplers,  filters 
and  transitions  [Ref.  41]. 
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The  geometric  configuration  of  the  step  discontinuity  and  its  equivalent  cir¬ 
cuit  are  shown  in  Figure  4.1. 
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Figure  4.1  Microstrip  step  discontinuity  and  the  equivalent  circuit 
The  effect  of  the  total  discontinuity  inductance  ( Ld )  may  be  separated  into  L\  and 
Lo  as 

U  =  ■  Ld  (H)  (4.1a) 

T  lJw2 

Lt  =  ,  LZ\  ^  (H)  (4.14) 

1  i  Ltw 2 

where  Lw\  and  Lw 2  are  the  inductances  per  unit  length  for  the  microstrip  of  width 
W\  and  W2  given  by 

Lwm  =  — m  =  1,2  (H/m)  (4.2) 

c 

where  Zom  is  the  characteristic  impedance  and  c«//m  is  effective  permittivity  of  the 
width  Wm,  and  c=3x  108  m/s.  The  closed-form  expressions  for  C<j  and  Ld  have 
been  derived  from  curve  fitting  numerical  results.  These  expressions  are  [Ref.  41] 

-^p=  =  (10.1  log  £,  +  2.33)[h  _  12.6  log  tr-  3.17  (pF/m) 

(for  er  <  10;  1.5  <  IFj/IFj  <  3.5)  (4.3a) 
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(nH/m)  (4.36) 


and 


where  tT  is  the  relative  dielectric  constant  and  h  is  the  substrate  thickness.  Equation 
(4.3a)  yields  a  percentage  of  error  less  than  10  percent  and  equation  (4.3b)  has  an 
error  less  than  5  percent  for  W1/W2  <  5  and  Wi/h  =  1.0. 

2.  Cascaded  Step  Discontinuity 

The  cascaded  microstrip  step  discontinuity  is  formed  by  combining  two  sin¬ 
gle  step  discontinuities  as  shown  in  Figure  4.2(a).  Short  (<  A?/4)  lengths  of  high 
impedance  (narrow  width)  microstrip  line  will  behave  predominantly  as  a  series  in¬ 
ductance  ( Ls ).  Whereas,  a  very  short  <C  Afl/ 4  length  of  low  impedance  (wide  width) 
line  will  act  predominantly  as  a  shunt  capacitance  (Ca)  [Ref.  18:  pp.  212-216]. 
End-inductances  (Le)  are  also  introduced  in  the  low  impedance  line  as  its  length 
approaches  a  quarter-wavelength.  The  predominant  series  inductive  reactance  of 
the  line  of  length  I2  is 

uLa  =  Z02  sinh  •  (4.4) 

Assuming  the  small-angle  approximation  such  that  (2nl2/Xg2  x/4),  equation  (4.4) 

is  approximated  as 


T  ^  ^02  h  , A 

L>  ~  Tv  (4-5) 

where  the  frequency-dependent  microstrip  wavelength  is 

A’2  =  /vA«//(/)' 

When  high  impedance  widths  ( W2)  are  equal,  the  equivalent  circuit  can  be  described 
by  Figure  4.2(b). 
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The  predominant  shunt  capacitive  suspectance  of  the  line  of  length 
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Figure  4.2  Cascaded  Step  Discontinuity  and  the  Equivalent  Circuit 
and  for  27t/i/Asi  <C  7r/4,  the  shunt  capacitance  is  approximated  by 

9  }  Zq\  Xgl 

Similarly,  the  end-inductances  are  obtained  as 

t  ~  Zoi  ll 
L'~ 


(4.7) 


(4.8) 
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The  relatively  low  characteristic  impedance  Z 01  will  cause  a  very  small  inductive 
effect  ( Lw\ )  at  each  step  discontinuity.  Therefore,  at  frequencies  up  to  a  few  giga¬ 
hertz,  the  inductance  L\  in  the  equivalent  circuit  can  be  neglected  [Ref.18:  p.  217]. 
A  modified  equivalent  circuit  is  formed  by  converting  the  low  impedance  equivalent 
T-network  into  a  ?r-network  as  shown  in  Figure  4.3(a).  The  shunt  capacitances 
Caj 2  and  the  discontinuity  capacitances  C<f  are  added  to  form  CT.  Also,  each  step 
discontinuity  inductance  Z-2  is  combined  with  their  adjacent  series  inductance  Ls  to 
form  a  new  inductance  T/,.  Assuming  very  short  microstrip  lengths  and  from  (2.24), 
the  new  inductance  can  be  approximated  as  a  lossless,  distributed  transmission  line 
section  having  a  characteristic  impedance  of  Z02  and  a  time  delay  approximated  by 


(4.9) 


Figures  4.3(b)  and  4.3(c)  show  the  complete  lumped-distributed  equivalent  circuits 
for  the  cascaded  microstrip  step  discontinuity. 


(a) 


(b) 

Figure  4.3  Modified  cascaded  step  discontinuity  equivalent  circuit 
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(c) 

Figure  4.3  Modified  cascaded  step  discontinuity  equivalent  circuit  (continued) 


3.  Capacitive  Shunt  Equivalent  Circuit  Model 

An  alternative  equivalent  circuit  is  presented  that  is  designed  to  yield  an  “all- 
pole”  transmission  network  function.  The  equivalent  ^-network  of  Figure  4.3(b)  is 
converted  into  a  T-network  and  the  series  inductances  are  combined  as  shown  in 
Figure  4.4(a).  The  combined  inductances  are  approximated  as  lossless,  distributed 
transmission  line  sections  having  time  delays  of 

(2  Lh  +  Le) 


t'd  = 


2Zo2 


(4.10) 


The  shunt  capacitance  equivalent  circuit  model  is  shown  in  Figure  4.4(b)  where  t'd 
represents  the  propagation  delay  time. 
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Figure  4.4  Capacitive  shunt  equivalent  circuit 
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Figure  4.4  Capacitive  shunt  equivalent  circuit  (continued) 

4.  Modeling  Propagation  Loss 

Chapter  II  presented  a  resistive  7r-network  attenuator  to  model  the  total 
propagation  loss  in  a  single  microstrip  line  section.  This  scheme  will  again  be  used 
for  the  equivalent  circuits  of  Figures  4.3(c)  and  4.4(b).  Figure  4.5  illustrates  how  the 
half-symmetrical  losses  of  the  cascaded  microstrip  step  discontinuity  are  modeled. 
Using  (2.8b)  and  (2.10),  a  total  propagation  loss  can  be  obtained  for  the  microatrip. 


Lossless  Reactive 
Equivalent  Network 


Figure  4.5  Lossy  cascaded  microstrip  step  equivalent  circuit 

Next,  we  equally  represent  the  total  loss  by  two  identical  resistive  7r- networks  as 
shown  in  Figure  4.5. 
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B.  PARAMETRIC  MODELS  FOR  MICROSTRIP  DISCONTINUITIES 
1.  Reflection  and  Transmission  Network  Functions 

The  characteristic  impedance  mismatch  caused  by  the  abrupt  change  in  the 
width  dimension  will  produce  both  reflected  and  transmitted  traveling  voltage  wave¬ 
forms.  The  objective  of  this  section  is  to  derive  reflection  and  transmission  network 
functions  which  accurately  characterize  the  cascaded  microstrip  step  discontinuity 
propagation  effects.  Initially,  a  third-order  transmission  network  function  is  derived 
from  the  equivalent  circuit  shown  in  Figure  4.3(c).  This  network  function  is  trans¬ 
formed  into  an  ARMA  model.  Next,  a  first-order  approximation  is  solved  from  the 
single  capacitive  shunt  equivalent  circuit  of  Figure  4.4(b).  This  single-impedance, 
single-capacitor  network  function  is  then  transformed  into  an  AR  model. 
a.  Third-Order  Transmission  Function 

The  reflection  network  function  for  the  7r-network  of  Figure  4.3(c),  ex¬ 
pressed  in  the  complex  s-domain,  is  defined  as  [Ref.  38:  p.21 7] 


r,  ,  loaM  -  YL(s) 

W  VtoW  +  YL(s) 


(4.11) 


where  the  admittances  >o2(-s)  and  Yi(s)  are  obtained  as  indicated  in  Figure  4.6. 


Figure  4.6  Characteristic  and  load  admittances 


37 


Substituting  for  Vo2(-s)  and  V/,(s)  in  (4.11)  and  manipulating  the  resultant  expression 
yields 


r  ,  .  = _ ICiL'Zfos3  +  (2 C,Zg2  H-  Lt)* _ 

'  1  (C2i,Z22)s3  +  (2C,L,Zm)s2  +  (2C,Zc22  +  L')s  +  2Zm 


(4.12) 


where  the  subscript  (3)  denotes  a  third-order  network  function.  The  transmission 
function  is  defined  as  [Ref.  3S:  p.  217] 


t3(s)  =  i  +  r  3(s). 


(4.13) 


Substituting  (4.12)  into  (4.13)  yields 


Tz(s)  = 


C,z02)  *2+  icjzl,)  S+(c2L'Zo2) 

2  V  |  (2C„2&  +  L'\  (  2 

C,Za2)  +V  C3L,Z22  )  \C;LtZm 


(4.14) 


b.  First-Order  Reflection  and  Transmission  Functions 

A  first-order  approximation  of  the  reflection  network  function  is  deter¬ 
mined  by  solving  (4.11)  for  the  equivalent  circuit  of  Figure  4.4(b).  Evaluating 
(4.11)  yields 


r  i(-s)  =  - 


s(CxZ02) 


\+s(CxZ02) 

and  the  transmission  network  function  becomes 


Ti(s)  = 


1  +  s  (CtZ02) 


(4.15a) 


(4.156) 


2.  ARMA/AR  Models  for  the  Cascaded  Step  Discontinuity 
a.  Autoregressive  Moving- Average  Parametric  Model 

The  first  step  is  to  perform  a  partial  fraction  expansion  of  equation  (4.14). 
The  s-domain  expansion  is  generalized  as 

,  toe-* 


T(  \  Kl  K2e  K2t 
3l5j  {S-dxy  {s-d2y  {s-d\) 


(4.16) 
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where  the  three  poles  are: 


dl  C*Z02 

d.2  =  —0  4-  ;w 

d*2  =  —0  —  jw- 


Using  the  impulse-invariant  design  method,  a  digital  filter  transfer  function  is  ob¬ 
tained  from  (4.16)  as 


TAz)  = 


I<  i 


1  -  e~dl  Tz~1 


+ 


2K2  cos (6)  -  2K2e  a  Tz  1  cos(fl  -  u  T) 
1  —  2e~a  T  cos (uj  T)z~l  +  e~d i  Tz~~. 


After  considerable  complex  algebra  the  generalized  ARMA  transfer  function  be¬ 
comes 


bo  +  b\z  1  -f-  b2Z  2 

1  —  Gl  2  — 1  —  G2  Z~2  —  azz~3 


(4.18) 


where 


bo 

bi 

b2 

ai 

«3 


K\  -f  2I<2  cos(0) 


-2 


I<2e  dlT  cos (0)  +  e  aT  {I<\  cos(u :T)  -f  I< 2  cos(0  -  uT)) 


- diT 


K\+2I\2£  aT  cos( 


s^-wT)] 

2e~aT  cos(uT)  -(-  e~dlT 

-  2e~{,T+dl)T cos{uT)  +  e~dlT 
—2d\T 


The  most  significant  filter  coefficient  03  identifies  the  time  constant  ( Zq2Ct )  of  the 
network.  In  terms  of  the  sampling  interval  T,  this  filter  coefficient  is 


03  =  exp 


( 


(4.19) 
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A  reflection  digital  transfer  function  is  obtained  from  the  bilinear  transformation  of 
equation  (4.15a)  as  [Ref.  34:  pp. 206-211] 


Ti(2)  =  ri(s)  2  /I  -  z~l 


=  2-(Lz£l\ 

TU  +  r>J 


Fi  (*)  =  - 


n  (z)  = 


2CtZ02(1  -z-1) 

(T  +  2CxZo2)  +  (r-2C,Zo2)--1 

1  +  a\z~l 


(4.20a) 


and  the  transmission  transfer  function  is  evaluated  as 


Ti{z)  = 


Ti(z)  = 


T(l  +  z -1) 

{T  +  2CxZ02)  +  (T  -  2CtZ02) 
6q(1  4-  z ~1) 

1  +  ajz-1 


(4.206) 


which  is  a  lst-order  ARMA  model.  The  parameters  are  defined  as  follows: 

bo==  ir+^CrZ 02  )’ 


Ll  _  (  2C*Zq2 
°~  \T  +  2CrZQ2)' 


=  (T-2CrZ02\ 
ai  \T  4-  2CxZ02/ 

6.  Autoregressive  Parametric  Model 

A  first-order  AR  digital  filter  of  the  transmission  network  function  is  ob¬ 
tained  using  the  matched  z-transform  design  method.  The  single  pole  of  equation 
(4.15b)  is  directly  related  to  the  pole  of  the  AR  transfer  function  [Ref.  42:  pp. 
666-667]  as 

(_Z_) 

r  _  \C.Zo2j  (,  01\ 


Tar(z)  = 


1  —  e  (c>r Z02}  Z-1 


(4.21) 
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c.  Realization  of  the  Lossy  Parametric  Model 

The  proposed  ARMA/AR  reflection  and  transmission  transfer  functions 
were  derived  from  lossless  reactive  networks.  The  propagation  losses  which  are  as¬ 
sociated  with  both  the  low  and  high  impedance  microstrip  lines  must  be  included  in 
the  parametric  models.  These  losses  are  represented  as  constant  “loss”  multipliers. 
The  value  of  the  multiplier  is  defined  by 

L  =  exp  [-(<*<  /)]  (4.22) 

where  /  represents  the  microstrip  line  length  and  a<  is  the  total  attenuation  factor 
in  nepers  per  unit  length.  Figure  4.7  describes  the  lossy  parametric  model  for  the 
cascaded  microstrip  step  discontinuity. 

3.  Parametric  Models  for  Multiple  Discontinuities 

In  this  section  a  parametric  model  is  developed  for  a  multi-section  microstrip 
structure  described  by  a  combination  of  cascaded  step  discontinuities.  Wave  prop¬ 
agation  in  the  multi-section  microstrip  structure  exhibits  several  of  the  physical 
characteristics  of  the  multi-layered  earth  model  developed  in  geophysical  signal 
processing  [Ref.  35:  pp.  306-315].  Similarly,  two  different  propagating  waves  exist 
within  each  microstrip  section,  an  incident  (right-propagating)  wave  and  a  reflected 
(left-propagating)  wave.  Unlike  the  lossless  multi-layered  earth  model,  propagating 
signals  in  microstrip  are  influenced  by  frequency-dependent  losses  and  dispersion. 
Energy  conservation  is  established  across  each  discontinuity  boundary.  However, 
energy  loss  will  occur  between  adjacent  discontinuities.  A  linear  time-invariant 
system  is  used  to  describe  each  discontinuity  boundary.  At  each  discontinuity  an 
incident  signal  is  partially  reflected  and  partially  transmitted  into  the  next  section. 
A  proportional  response  of  each  component  is  given  by  the  concurrent  convolution 
of  the  incident  signal  with  both  a  reflected  and  a  transmitted  impulse  response  of 
the  linear  system.  A  lattice-like  system  flowgraph  relating  right-propagating  t^(r), 
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and  left-propagating  vm(z)  signals  at  the  m-th  discontinuity  boundary  is  shown  in 
Figure  4.8.  The  reflection  transfer  function  rm(z) 


Figure  4.7  lst-order  ARMA  model  realization  of  the  cascaded  step  discontinuity 

and  the  transmission  transfer  function  Tm(z)  are  the  lst-order  ARMA  models  of 
(4.20a)  and  (4.20b),  respectively.  The  overall  parametric  model  for  multiple  dis¬ 
continuities  may  then  be  constructed  by  using  the  lattice  as  building  blocks,  and 
interconnecting  the  lattice  sections  with  lossy  delay  lines.  Figure  4.9  describes  two 
adjacent  discontinuity  sections  of  a  multi-section  structure.  Here  an  incident  voltage 
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sequence  t4n"^[n]  initially  enters  a  lossy  delay  line  Dm{z),  and  a  reflected  voltage 
sequence  dL  ^[n]  is  the  resultant  outcome  after  multiple  reflections  from  the  other 
discontinuities. 
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Figure  4.8  Discontinuity  lattice  configuration 

The  flowgraph  of  the  model  shows  an  interaction  of  right-  and  left-propagating 
waves  through  the  multi-section  structure.  We  will  characterize  the  lossy  delay  line 
for  an  m-th  section  by  the  transfer  function 

Dm{z)  =  e~tQ<  }  *"»)  z~Tm  (4.23) 

where  a[m ^  is  the  m-th  section  attenuation  factor,  lm  the  section  length,  and  Tm  is 
the  discrete-time  propagation  delay  per  section. 

In  summary,  both  equivalent  circuit  and  ARMA/AR  parametric  models  were 
presented  for  a  cascaded  microstrip  step  discontinuity.  Losses  were  included  in 
the  circuit  model  by  resistive  7r-networks.  Corresponding  lossless  reflection  and 
transmission  network  functions  were  derived  from  the  reactive  networks.  Next, 
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discrete  transformations  were  performed  on  the  network  functions  to  derive  the 
ARMA/AR  digital  transfer  functions.  Multiple  microstrip  step  discontinuities  were 
modeled  by  cascading  the  lossy  equivalent  circuit  of  a  single  step  discontinuity. 
A  lattice-like  lst-order  ARMA  model  was  introduced  to  describe  a  multi-section 
microstrip  structure.  In  Chapter  VI,  transient  analysis  simulations  of  the  equivalent 
circuit  models  will  be  compared  with  experimental  measurements  to  validate 


Figure  4.9  Multi-section  lattice  model 

the  accuracy  of  the  circuit  models.  In  the  next  chapter,  we  will  develop  several  pa¬ 
rameter  estimation  algorithms,  and  introduce  new  system  identification  techniques 
for  multi-section  microstrip  structures. 
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V.  PARAMETER  ESTIMATION  ALGORITHMS 


The  parameter  estimation  algorithms  presented  are  based  upon  deterministic 
and  stochastic  methods.  Figure  5.1  describes  an  unknown  system  to  be  identified. 
The  observable  input  and  output  signals  can  be  either  deterministic  or  stochastic. 
The  ARMA  and/or  AR  model  parameters  of  the  unknown  system  are  determined  by 
the  parameter  estimation  algorithms.  Additionally,  simulation  results  are  presented 
to  study  the  performance  of  the  algorithms. 


Figure  5.1  Block  diagram  representation  of  the  p'  uneter  estimation  problem 

The  deterministic  methods  to  be  presented  are  the  weighted  least  squares  (WLS) 
algorithm,  the  network  function  approximation  method,  and  the  Schur  algorithm. 
The  WLS  algorithm  is  a  method  for  finding  the  filter  coefficients  of  a  p-th  order 
recursive  digital  filter,  which  gives  an  optimum  least  squares  approximation  to  a 
known  sample  impulse  response.  The  network  function  method  approximates  an 
s-domain  network  function  from  the  unit  step  response  characteristics,  such  as  rise 
and  delay  times,  of  a  propagating  pulse.  Finally,  the  Schur  algorithm  estimates 
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the  AR  model  parameters  from  a  given  sample  impulse  response.  Our  emphasis 
will  focus  on  the  deterministic  methods  because  avaliable  microwave  measurement 
techniques  have  deterministic  excitation  sources. 

In  the  discussion  of  the  stochastic  methods  we  assume  that  a  nearly  white 
wideband  microwave  noise  source  excites  the  microstrip  structure.  When  the  noise 
source  and  the  output  random  data  are  both  available  for  signal  processing,  an 
ARMA  parameter  estimation  algorithm  based  on  a  generalized  Mullis- Roberts  (M- 
R)  criterion  is  used  [Ref.  15].  However,  if  the  input  noise  source  cannot  be  measured, 
an  estimate  of  the  driving  source  can  be  obtained  from  the  output  process.  This 
restriction  in  observable  data  precludes  the  direct  application  of  the  general  estima¬ 
tors.  Therefore,  a  suboptimal  two-stage  least  squares  ARMA  estimation  algorithm 
will  be  presented  to  solve  this  problem. 

In  Chapter  IV,  we  presented  a  lattice  parametric  model  for  the  multi-section 
microstrip  discontinuity.  New  algorithms  will  now  be  introduced  that  will  estimate 
the  ARMA/AR  parameters  for  this  model  using  reflected  and  transmitted  time- 
domain  measurements. 

A.  DETERMINISTIC  PARAMETER  ESTIMATION 

1.  Weigthed  Least  Squares  (WLS)  Algorithm 

a.  WLS  Algorithm  Derivation 

Evans  and  Fischl  [Ref.  43]  presented  a  computational  algorithm  for  es¬ 
timating  the  coefficients  of  a  pth-order  ARMA  digital  filter  whose  unit  impulse 
response  best  approximates,  in  a  minimum  mean-square  error  sense,  a  prescribed 
impulse  response  of  finite  length.  Here  we  derive  a  slightly  modified  version  of  their 
original  least  squares  solution  to  estimate  the  ARMA  parameters. 

The  transfer  function  Harma(z)  of  (3-2)  can  be  expressed  in  terms  of  samples 
of  the  impulse  response 

HaRMa{z)  =  ho  +  h\z~l  +  h2z~2  +  • . .  •  (5.1) 
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The  first  m  samples  of  the  estimated  impulse  response  forms  the  m-vector 

ho 

A  -  hi 


and  the  m  samples  of  the  prescribed  impulse  response  is 


The  error  vector  is  defined  by 


e  =  h  —  h 


By  long  division  of  (3.2),  we  can  relate  the  AR  and  MA  parameters  (a,  6)  to  the 
estimated  impulse  response  h(a,b )  as 


(5.5a) 


where 


and  a  — 


(5.5  b) 


m  = 


ho  0 

hi  ho 

A  *  A  * 

.ftp—]  ftp— 2 


m  = 


hp  hp-i 
_  hp+ j  hp 


htn— 1  hm—2 


0  0 
0  0 

ho  0 


hm—p—  1 


(5.5c) 


(5.5d) 


A  least  squares  solution  can  be  determined  by  minimizing 


II  e  (<z,  6)  ||  =  ^e,(a,6) 
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with  respect  to  the  coefficients  (a,  6). 

Burrus  and  Parks  [Ref.  44]  have  proposed  a  linear-equation  error  vector  in  the 
matrix  formulation,  given  by 

d{a)  =  [H2]  a  (5.7 a) 


where 


d(a)  = 

I 

djn—p{a ) 

The  error  vector  e(a,b)  of  (5.4)  is  related  to  d(a)  by  rearranging  (5.7a)  as 


d\(a) 

d2(a) 


(5-7  b) 


d(a)  =  [ B{a)]T  h 

where  the  m  x  (m  —  p)  matrix 


Qp 

0 

...  0  1 

O.p—1 

dp 

...  0 

B(a)  = 

a  o 

a\ 

. . .  dp 

0 

do 

...  d\ 

.  0 

6 

...  a0. 

(5.8) 


(5.9) 


We  can  substitute  h  of  (5.4)  into  (5.8),  and  recognizing  from  the  lower  partition  of 
(5.5a)  that 

[B{a)]T  h  =  0,  (5.10) 

we  have 

d(a)  =  [B(a)]T  e(a,  b).  (5.11) 


An  inverse  relationship  between  the  errors  is  formed  using  a  transformation  weight¬ 
ing  matrix  [W(a)\.  Hence,  the  minimization  problem  of  (5.6)  can  now  be  expressed 
as  [Ref.  44] 

min  j|  e(a,  b)  ||  =  min  ||  [W{a)]d{a)  ||  (5.12) 

a,  b  a 

where 

[H^(a)]  =  [H(a)][H(a)rH(a)]_1.  (5.13) 
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The  computational  algorithm  begins  with  an  initial  estimate  of  Equation 
(5.7a)  can  be  partitioned  as 


d(a)  =  [-/  |  [X]]  a 


where 


m  = 


hp— i 

lip  .  .  , 


ho 

hi 


[_  ^771  —  2  *  *  *  hm—p—\ 

hn 


f  = 


and  the  AR  parameter  vector 


hp+ 1 

hm—l 


a  i 
<*2 


a  = 


The  least  squares  solution  of  a  yields 

a(0)  =  [xT  x]~l  xT  f 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


The  initial  estimate  equally  weights  the  error  vector  d(a).  However,  the  error 
formulation  can  be  weighted  differently  depending  upon  the  initial  estimate.  Then, 
the  weighted  least  squares  solution  at  the  z-th  iteration  is  given  by 

XTW  x\ 


fiW  = 


XTW  / 


(5.19) 
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where  the  weighting  matrix  j^VU  j  is  in  terms  of  the  previous  estimate  a)1-1). 

The  iteration  is  continued  until  a ^  Rs  Using  the  optimal  estimate  for  a,  the 

weighted  error  eu,(a,6)  is  calculated  by 

s.M)=  [w  (a(,))]  d  (a«) .  (5.20) 

From  (5.4),  the  estimate  h(a,b)  is  given  by 

h(a,b)  =  h  —  ew(a,b).  (5.21) 

Finally,  the  numerator  coefficient  estimate  vector  b  is  calculated  directly  by 

b=[Hl]a.  (5.22) 


b.  Simulation  Results 

The  weighted  least  squares  algorithm  was  implemented  using  a  Fortran 
program.  Several  reference  models  were  estimated  beginning  with  an  ARMA  order 
(3,2),  denoting  AR  order  p  =  3,  and  MA  order  q  =  2.  The  transfer  function  is 


H(z) 


0.5  —  0.4z-1  +0.892-2 
1  -  0.202-1  -  0.25z-2  +  0.05;~3 ' 


(5.23) 


The  actual  reference  model  parameters  and  the  WLS  algorithm  parameter  estimates 
are  shown  in  Table  1.  All  results  are  obtained  after  one  iteration  of  the  algorithm 
using  20  samples  of  the  impulse  response. 

Table  1.  WLS  SIMULATION  (3,2)  MODEL  RESULTS 


Parameter 

Actual 

Estimated 

AR: 

0.2000 

0.2500 

-0.0500 

-0.0500 

MA: 

0.5000 

-0.4000 

0.8900 

0.8899 
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We  next  consider  a  second  reference  model  with  order  (4,2)  having  transfer  function 


H{z) 


1  +  0.2z-1  -0.992-2 

1  -  0.20.J-1  +  0.62^-2  -  0.152c-3  +  0.3016c-4 


and  the  obtained  results  are  given  in  Table  2. 


(5.24) 


Table  2.  WLS  SIMULATION  (4,2)  MODEL  RESULTS 


Parameter 

Actual 

Estimated 

AR: 

0.2000 

0.2002 

-0.6200 

-0.6197 

0.1520 

0.1529 

-0.3016 

-0.3014 

MA: 

1.0000 

1.0000 

0.2000 

0.1998 

-0.9900 

-0.9902 

The  final  example  is  a  strictly  AR  model  of  order  (4,0)  having  a  transfer  function 


H(z) 


1 

1  -  0.20r-1  +  0.6222'2  -  0.1512-3  +  0.3552~4  ' 


(5.25) 


The  actual  and  estimated  AR  parameters  are  listed  in  Table  3. 

Table  3.  WLS  SIMULATION  (4,0)  MODEL  RESULTS 


Parameter 

Actual 

Estimated 

AR: 

0.2000 

0.2004 

-0.6220 

-0.6221 

0.1510 

0.1513 

-0.3550 

-0.3551 

The  above  examples  demonstrate  the  estimation  accuracy  of  the  WLS  algorithm 
for  (q  <  p)  order  models. 

2.  Network  Function  Approximation  Method 

The  following  method,  introduced  by  Elmore  [Ref.  45],  can  approximate 
the  normalized  network  function  from  the  unit  step  response  delay  and  rise  times. 
Assume  that  the  microstrip  section  is  excited  by  a  unit  step  and  that  its  step 
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response  v2(t),  has  been  normalized  as  t>2n (t).  The  microstrip  group  delay  time  is 
then  defined  as 

OC 


J  t  *4 n(t)dt 

- - =  J  t  hn(t)dt  (5.26) 

/  v'2 n{t)dt  0 


where  v'2n (t)  =  hn(t)  is  the  normalized  impulse  response.  The  rise  time  is  defined 
as 


Tr  =  V2n 


=  V2n 


OO 


J  (t  -  Tg)2  hn(t)dt 


1/2 


Tt 


OO  OO  OC 

J  t2  hn(t)dt  -2  Tg  J  t  hn(t)di  +  t2  J  hn[t)dt 
.0  0  0 
OO 

J  t 2  hn 


{ t)dt  -  t2 


0 

1/2 


1/2 


(5.27) 


Figure  5.2  describes  a  network  voltage  output  response  resulting  from  a  unit  step 
input.  The  normalized  network  function  is  expressed  as 

OO 

Hn{s)  =  J  hn(t)e~sidt  (5.28) 

o 


Expanding  e  J  in  a  power  series  in  (5.28)  yields 

on 

„2  j2 


#"(•$)  =  J  hn{t){\  -St  +  - .  .  .)dt, 

0 

oo  oo  2  00 

=  J  hn(t)dt  -  s  J  t  hn{t)dt  +  ^  J  t2  hn(t)dt  - - 

oo  o 

Kg-?)-- 


=  1  -  ST‘  +  2! 


(5.29) 
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As  shown  in  Figure  5.3,  the  microstrip  transmission  line  is  described  by 
M  identical  network  sections  in  cascade  each  with  independent  propagation  delay 


Figure  5.2  A  network  N  with  its  unit  step  response 

times  Tgi ,  rg2 , . . . ,  TgM ,  and  rise  times  Tri ,  Tn ,  . . . ,  Tru .  Then  the  overall  normalized 
network  function  is 
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Finally. 


M  2  f  M  t2  (  M  \  ' 

tfn(s)  =  i-s£v  +  |r  Z^r+(ZT^)  •••  (5-31^ 

Jb=l  '  U= 1  \Jb=l  /  . 

where  the  total  propagation  group  delay  is 

M 

r,  =  Y,  V  (5-32) 

fc-1 

and  the  overall  rise  time  for  M  sections  is 

M 

7V2  =  £r,V  (5.33) 

k=  1 


Figure  5.3  Microstrip  transmission  line  as  M  cascaded  networks 

Assume  that  the  propagating  pulse  is  the  response  of  some  networ1'  to  an  unit  step. 
Evidently,  the  rise  time  of  the  previous  network(s)  is  7Vi .  Now  the  pulse  continues 
through  an  unknown  network  which  has  a  rise  time  Tr.  The  resulting  pulse  then  has 
a  total  rise  time  Tr 2-  FYom  (5.33),  the  expected  rise  time  of  the  unknown  network 
is  given  by 

Tr  =  Jt?2  -  7?,.  (5.34) 
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After  estimating  the  network  delay  and  rise  times,  the  ARMA  filter  coeffi¬ 
cients  of  equation  (3.11)  are  solved  from  the  polynomial  coefficients  of  (5.29),  while 
the  first-order  AR  filter  coefficient  is  also  related  to  rg  and  Tr  as  follows: 


a  =  exp 


where  T  denotes  the  sampling  interval. 


(5.35) 


3.  Deterministic  Schur  Algorithm 

The  autocorrelation  function  (ACF)  of  the  normalized  impulse  response  of 
equation  (3.16)  is  computed  as 


,°° 

R(t)  =  I  h(t)h(t  +  r)dt,  for  r  =  0, 1, 2, ...,  N  (5.36) 

Jo 

where  r  is  the  correlation  lag.  However,  the  sampled  autocorrelation  function  will 
give  a  discrete  correlation  vector  as: 


[R(0),R(1),R(2),  ...  ,R(N)]  (5.37) 


Next,  the  ACF  data  is  applied  to  the  Schur  algorithm.  The  Schur  algorithm  is  a 
recursive  algorithm  that  solves  the  AR  model  parameters  (also  called  the  partial 
reflection  coefficients)  [Ref.  14].  The  algorithm  begins  by  forming  a  generator 
matrix  ( Gq  )  using  the  sampled  autocorrelation  vector  as 


\R(0) 

...  R(N) 

0 

R(l) 

1  R{  2) 

...  R(N) 

Shifting  the  first  column  down  yiel  is 


0 

R(0) 

R(l) 

0 

R(l) 

m 

R(N-  1) 
R(N) 


(5.38) 


(5.39) 
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The  lst-order  model  reflection  coefficient  is  computed  as  the  ratio  of  the  (2,2)  and 
(2,1)  terms  of  (5.39)  as 

kl  =  I(of  (5'4°) 

A  new  generator  matrix  is  then  formed,  given  by 


cf  =  e(t,  )Gj 


(5.41) 


where 


Similarly,  the  2nd-order  model  reflection  coefficient  is 

R(2)~k1R(l) 

2  R(0)~kiR(iy 


(5.42) 


(5.43) 


The  algorithm  repeats  until  the  desired  order  reflection  coefficient  is  obtained. 
When  the  impulse  response,  h(t),  describes  a  lst-order  AR  transfer  function,  a  sin¬ 
gle  reflection  coefficient  will  result.  From  this  reflection  coefficient  the  time  constant 
of  the  impulse  response  is  easily  found.  The  computational  efficiency  is  achieved  by 
iterating  the  autocorrelation  vector  through  an  Nth-order  AR  lattice  structure.  In 
the  next  section,  the  Schur  algorithm  will  be  applied  to  an  AR  stochastic  process. 


B.  STOCHASTIC  PARAMETER  ESTIMATION 

The  following  algorithms  are  presented  to  estimate  the  ARMA/AR  model  pa¬ 
rameters  when  the  observable  input /output  signals  are  stochastic  processes.  Assum¬ 
ing  a  white  noise  process  as  the  driving  input,  the  system  to  be  identified  produces 
an  output  process  as  shown  in  Figure  5.4.  The  efficiency  of  a  particular  estimation 
algorithm  is  governed  by  the  choice  of  model  selected  and  its  order.  In  practice,  we 
do  not  usually  know  a  priori  which  model  to  choose.  Once  a  model,  either  ARMA 
or  AR,  has  been  chosen,  we  must  specify  the  model  order.  In  choosing  a  model  for 
the  microstrip  transmission  line,  we  will  select  the  ARMA/AR  models  proposed  in 
Chapter  III. 
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Autoregressive  models  are  the  most  widely  used  models  because  the  analysis 
algorithms  for  extracting  the  model  parameters  are  found  by  solving  a  set  of  linear 


Figure  5.4  Stochastic  system  identification  problem 


equations.  When  the  AR  modeling  assumption  is  valid,  these  algorithms  provide 
very  good  estimates  of  the  model  parameters  [Ref.  33:  p.  131].  Unfortunately, 
application  of  these  algorithms  to  non-AR  time  series  data  usually  results  in  poor 
quality  estimates.  Furthermore,  if  an  AR  estimator  is  applied  to  a  process  that  is 
not  AR,  then  the  true  AR  model  would  be  one  of  infinite  order.  Any  finite  order 
AR  model  will  introduce  bias  errors  from  modeling  inaccuracies.  A  trade-off  must 
take  place  to  choose  order-p  large  enough  to  reduce  the  bias  or  to  choose  order-p 
small  enough  to  reduce  the  estimation  errors. 

In  ARMA  modeling,  the  best  least  squares  estimate  of  the  model  parameters 
is,  generally,  a  nonlinear  function  of  the  past  observations.  Nonlinear  optimiza¬ 
tion  techniques  are  usually  computationally  intensive  and  may  not  converge  to  the 
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global  minimum.  Two  ARMA  parameter  estimation  algorithms  will  be  presented 
for  stochastic  data. 

1.  Schur  Algorithm  -  Revisited 

The  AR  filter  coefficients  (ai, 02,03,  ...,ap)  are  related  to  the  autocorrela¬ 
tion  matrix  of  the  data  to  be  modeled  and  the  noise  variance  (a2)  by  the  normal 
equations: 

Ryy(0)  Ryy(-l)  .  Ryy{—n  +  l) 

Ryy(l)  Ryy(O)  .  Ryy(n) 

_Ryy{n—\)  ...  Ryy(l)  Ryy(0) 

Several  algorithms  have  provided  efficient  methods  for  solving  the  normal  equations 
[Ref.  33:  pp.  116-118].  The  Levinson  algorithm  for  solving  Toeplitz  systems  of  lin¬ 
ear  equations  is  well-known  in  digital  signal  processing.  However,  the  Levinson 
algorithm  has  been  discovered  to  be  computationally  less  efficient  than  an  alter¬ 
native,  the  so-called  Schur  algorithm,  because  of  a  requirement  to  evaluate  vector 
inner  products  [Ref.  12:  pp.  6-22]. 

We  have  previously  introduced  the  Schur  algorithm  for  estimating  AR  model  pa¬ 
rameters  using  deterministic  signals.  Assuming  that  the  unknown  system  is  driven 
by  a  white  noise  process,  the  AR  model  parameters  can  be  computed  from  the 
output  process  using  the  Schur  algorithm.  The  time  average  approximation  of  an 
autocorrelation  function  (called  the  sample  autocorrelation)  is  defined  as  [Ref.  46: 
p.  27] 

N-l-k 

Rvy[n]  =  jy  ^2  for  0  ^  n  ^  (A'  ~  !)  (5.45) 

Jt=o 

where  t/[n]  are  the  measured  output  signal  samples  of  length  N.  A  derivation  of  the 
Schur  algorithm  for  stochastic  signals  is  presented  in  Appendix  A. 

2.  Generalized  Mullis-Roberts  (M-R)  Algorithm 

An  ARMA  parameter  estimation  algorithm  based  on  a  generalized  Mullis- 
Roberts  criterion  has  been  proposed  by  Miyanaga  et  al.  [Ref.  47:  pp.  619-621], 


1 

ro 

a\ 

= 

*  O 

.  R, 

_ 1 

- 1 

..  O 

_ 1 

(5.44) 
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This  algorithm  is  derived  by  forming  forward  and  backward  prediction  models  using 
available  input  noise  and  output  data.  Prediction  errors  are  formed  to  satisfy  specific 
orthogonality  conditions  which  are  similar  to  those  found  in  AR  modeling  problems 
[Ref.  47:  pp.  116-121].  Finally,  autocorrelation  and  crosscorrelation  functions 
are  calculated  and  then  the  forward  and  backward  ARM  A  filter  coefficients  are 
predicted  from  the  correlation  lags.  A  previous  reference  model  is  estimated  using 
the  M-R  algorithm  in  order  to  compare  parameter  estimates  obtained  by  the  WLS 
algorithm.  The  transfer  function  of  order  (4,2),  given  by  (5.24),  was  driven  by  a 
white  noise  sequence.  Table  4  shows  the  ARMA  parameter  estimates  obtained 


Table  4.  M-R  SIMULATION  (4,2)  MODEL  RESULTS 


Parameter 

Actual 

(WLS)  Estimate 

(M-R)  Estimate 

AR: 

0.2002 

0.1933 

-0.6197 

-0.6165 

1 

0.1529 

0.1462 

-0.3014 

-0.3010 

MA: 

1.0000 

| 

0.1998 

I 

-0.9900 

-0.9902 

-0.9893 

from  the  M-R  algorithm. 

3.  A  Two-Stage  Least  Squares  Algorithm 

When  the  input  noise  is  unobservable,  the  M-R  algorithm  cannot  be  applied 
directly.  A  modified  two-stage  least  squares  approach  is  now  presented.  The  name 
two-stage  least  squares  is  due  to  the  use  of  a  large  order  AR  least  squares  estima¬ 
tor  followed  by  another  ARMA  parameter  estimation  algorithm.  This  approach 
is  known  to  be  suboptimal  because  the  noise  input  is  unobservable  and  must  be 
estimated  [Ref.  33:  p.  320].  Nevertheless,  we  modify  the  two-stage  least  squares 
algorithm  to  estimate  the  specific  ARMA  model  parameters  of  equation  (3.11).  The 
modification  uses  the  M-R  algorithm  instead  of  the  Levenison  algorithm  as  sum¬ 
marized  in  Table  6.  A  reference  model  was  estimated  having  a  transfer  function 
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of 


„  0.11722"1  . 

HaRMA  =  l-1.45i3,-iT0.57bF5  (5'46> 

at  a  sample  period  of  1.25  ps.  Table  5  lists  the  results  of  this  algorithm.  The  M-R 
algorithm  (Step  3)  produced  very  poor  MA  order  coefficient  results  from  the  noise 
estimate.  However,  the  AR  order  coefficients  ai  and  ct2  were  used  to  calculate  the 

Table  5.  TWO-STAGE  ALGORITHM  RESULTS 


Parameter 

Actual 

Estimated 

bi 

0.1172 

0.1097 

a  i 

1.4513 

1.4591 

Q2 

-0.5700 

-0.5701 

single  MA  order  coefficient  with  favorable  results. 

C.  MULTI-SECTION  MODEL  PARAMETER  ESTIMATION 

In  this  section  new  algorithms  are  presented  for  the  identification  of  the  multi¬ 
section  microstrip  structure.  The  layer  reflection  coefficient  estimation  algorithm 
will  approximate  the  number  and  location  of  cascaded  microstrip  step  discontinuties 
from  the  measured  reflection  data.  The  layer-probing  algorithm  is  an  iterative 
impulse  response  estimation  technique.  The  algorithm  assumes  that  the  multi¬ 
section  microstrip  structure  can  be  described  by  the  lattice  AR  parametric  model 
of  Chapter  IV. 

1.  Layer  Reflection  Coefficient  Estimation  Algorithm 

A  single  discontinuity  between  two  lossy  transmission  lines  having  different 
values  for  the  characteristic  impedance  Zm,  and  propagation  function  7m  is  shown  in 
Figure  5.5.  The  m-th  subscript  denotes  a  single  layer  defined  by  equal  propagation 
time  delay.  Assuming  that  a  current  wave  iffl'  is  incident  at  the  discontinuity  from 
the  m-th  layer,  we  find  that  the  incident  current  wave  alone  cannot  satisfy  the 
boundary  conditions  at  the  junction  [Ref.  35]. 
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Table  6.  SUMMARY  OF  TWO-STAGE  LEAST  SQUARES  ALGORITHM 

Given  the  output  process,  y[n],  of  length  A,  the  sample  interval,  T, 
and  L  N 

Step  1:  Predict  the  Tth-order  AR  parameters  ( a *) 

Using  the  Levinson  algorithm  [Ref.  33:  pp.  213-214]  with  output,  (y[n]) 

Step  2:  Estimate  the  noise  input  ( v  M) 

L 

u[n]  = 

fc=o 

Step  3:  Estimate  AR  filter  coefficients  (ai)  and  (02) 

Use  the  Mullis-Roberts  criterion  algorithm  [Ref.  15] 
where  (u[?i])  =»  input  and  ( y[n ])  output 

Step  4:  Solve  ((u;0  T) 

Substitute  (ai)  of  Step  3  =£•  (3.11c) 

Step  5:  Solve  (u>o\/l  —  C2  T) 

Substitute  (ai)  from  Step  3  =>  (3.11b) 

Substitute  (£wo  T)  from  Step  4  =>  (3.11b) 

Step  6:  Solve  MA  filter  coefficient  (6j) 

Use  ((wo  T)  from  Step  4,  and 
(uoy/l  —  C2  T)  from  Step  5  to  solve  £ 


Hence,  a  reflected  current  wave  iL  ^  is  generated.  Using  the  boundary 
at  the  impedance  mismatch,  we  then  write 


W + .t0'  =  4+i 


+  i<-> 

'  *m+l 


where 


Figure  5.5  Impedance  discontinuity  equivalent  circuit 
The  voltage  and  current  wave  relationships  are  defined  as  follows: 


.(+)/  _  ti+)' 

lm  — 

6m 

(_),  _  vfc0' 

lm  —  “^7 

6m 

,,(+) 

■(+)  _  vm+l 

*m+l  ~ 


j(~)  = 

‘m+1 


Zm+ 1 
(-) 


V, 


m-f  1 


Zm+1 


conditions 

(5.47a) 

(5.47fc) 

(5.47c) 


(5.48a) 

(5.48b) 

(5.48c) 

(5.48d) 
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(5.49) 


By  substituting  (5.48)  into  (5.47)  we  have 

m  m  _zm+1rra+!  w+v- 

The  instantaneous  voltage  at  the  mismatch  is  defined  as 


-  .,(+)'  _ 

m 


v(t,lm)  =  v\n>  -V, 

_.,(+)  _J-) 


Taking  (5.49)  and  (5.50),  we  solve  for  the  incident  voltage  wave  as 

-) 

JmJ 


„(+>'  -  J+> 

hm  ~  um+l 


+  ( En 

+  lm+l  1 

\  Tm 


where  the  reflection  coefficient  pm,  and  transmission  coefficient  t„ 
are 


Pm  — 


and 


Tm  — 


Zm+l  Zr 
Zm+1  4"  Zn 

2  Zm 


Zm  + 1  +  Zm 

Tm  ~  1  4"  pm- 

While  the  reflected  voltage  wave  at  the  m-th  layer  is 

(_),  (+)  (Pm\  (-)  f  J_ 

m  "  m+1UJ  m+1U» 

=  vm  ^  (Pm)  ~  vm+ 1  0  —  Pm)- 

Using  (5.51a),  the  voltage  wave  incident  at  the  ( m  +  l)-th  layer  is 

(!  +  Pm)  ~  Vmh  (pm)- 


Combining  (5.52)  and  (5.53)  into  a  scattering  matrix  yields 


r^(+) " 

1  4”  Pm 

Pm 

UrN 

Pm 

1  Pm 

'  „(+)'  ' 

Vm 

J-) 

_  ^  m  +  ] 

(5.50) 


(5.51a) 
at  the  mismatch 

(5.516) 


(5.51c) 


(5.52) 


(5.53) 


(5.54) 
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Next,  we  normalize  the  voltage  waves  with  respect  to  their  characteristic  impedance 


as 


and 


fvL+),l 

l 

v(-)/ 

_  vm 

V%n 

(5.55a) 


v"-6 1 

i 

|V+)  1 

vm+ 1 

.vL+i. 

VZm+l 

v(~) 
L^ro+l  J 

(5.556) 


Using  (5.54)  and  (5.55),  the  normalized  forward  recursion  matrix  yields 


rv<+>  1 

vm+l 

1 

'  1 

~  Pm 

1 

< 

3± 

-« 

v<-> 
vm+l  J 

V1  ~  Pm 

Pm 

1 

lT»j 

(5.56) 


The  normalized  incident  and  reflected  voltages  in  (5.56)  can  be  related  to  the  prop¬ 
agation  function  (7m)  as 


v 
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(+>' 

m 

(-)' 

m 


e  lm  lm  0 
0  elm  lm 


V 

V 


(+) 

m 

(-) 


g  —  (  Ck  m  0m  )  l  m  Q 

0  €(°m+]0  m)lm 


1 

’„(+)' 

Vm 

,  (-) 

Vm 

(5.57) 


Therefore,  (5.56)  can  be  expanded  as 


Vm+1 

1 

1  Pm 

e~°m  lm  0 

.  vm+l  . 

\/i  -pI 

Pm  1 

0  e°m  ‘m 
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m  *m  Q 
0m  Jm 


V 
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(+) 


(5.58) 


For  a  lossy  dispersive  microstrip  transmission  line  the  attenuation  am.  and  phase 
l3m  terms  are  functions  of  frequency. 
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Hence,  the  total  attenuation  per  m-th  layer  is  described  by 


Or n(f)  =  <*c(f)  +  Qrf(Z) 


(5.59) 


where  the  conductor  loss  ac,  and  the  dielectric  loss  aj  were  given  by  (2.8b)  and 
(2.10),  respectively.  Whereas,  the  phase  term  is  related  to  the  frequency-dependent 
effective  microstrip  permittivity  e£//(/)  as 


M)  = 


(5.60) 


Both  the  attenuation  and  phase  terms  are  approximated  as  constants  at  the  m-th 
layer.  Furthermore,  the  phase  term  0m  is  proportional  to  the  m-th  layer  propagation 
delay  time  [Ref.  35]  as 


t-m  — 


J9m 


=  /« 


d&m 

duj 


(5.61) 


where  vgm  is  the  group  velocity  in  the  m-th  layer. 

Next  we  develop  the  discrete-time  recursion  equations  that  solve  the  inverse 
scattering  problem.  Talcing  the  z-transform  of  the  recursion  matrix  (5.58)  and  using 
(5.61)  yields 


vL+i(z)' 
yj  i(z). 


VL+)(Z) 

.vL_)(z). 


(5.62) 


where 


,-i 


0  -  Pm 


Z  1  -  Pm 

PmZ  ^  1 


—  Q'm  In 


The  variable  z~l  represents  the  two-way  propagation  time  delay  of  the  m-th  section, 
while  z~ 5  represents  the  one-way  time  delay.  Using  (5.58),  the  (m  +  l)-th  layer 
reflection  coefficient  is  computed  by 


Pm  + 1 


v(-\ 
v  m-H 

0T 

vm  +  l 


(5.63) 
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A  layer  characteristic  impedance  profile  (in  ohms)  is  determined  from  the  layer 
reflection  coefficient  sequence  as 


'm+ 1 


= 


1  A  Pm  + 1 
.  1  Pm+1  . 


(5.64) 


where  the  initial  impedance  Zq  is  the  measured  input  characteristic  impedance.  In 
order  to  compensate  for  losses,  the  algorithm  must  have  prior  knowledge  of  the 
total  attenuation  per  layer  am.  Unfortunately,  this  loss  information  is  not  available 
for  the  inverse  scattering  problem,  but  the  overall  attenuation  per  structure  can  be 
estimated  from  the  measured  incident,  reflected,  and  transmitted  signals.  Under 
this  assumption,  we  can  obtain  an  average  attenuation  factor  for  each  layer.  A 
method  to  estimate  the  attenuation  factor  from  measured  results  will  be  introduced 
in  the  last  section. 

2.  Layer-Probing  Algorithm 

The  multi-section  microstrip  structure  is  modeled  as  cascaded  linear  time- 
invariant  systems.  Figure  5.6  describes  the  first  iteration  of  the  layer-probing  algo¬ 
rithm.  The  algorithm  begins  by  exciting  the  structure  with  an  impulsive  signal  vs[n). 
A  minimum-phase  transfer  function  T\{z)  is  estimated  for  the  first  discontinuity  us¬ 
ing  the  WLS  algorithm.  Next,  an  inverse  filter  A\(z)  of  the  transfer  function  T\(z)  is 
implemented  and  excited  by  the  original  impulsive  signal.  The  inverse  filter  output 
now  becomes  a  new  excitation  signal.  This  new  input  is  convolved  with  the  first 
discontinuity  transmission  impulse  response,  resulting  in  an  impulsive  signal  t>p[n] 
which  probes  the  second  discontinuity.  The  subsequent  reflected  wave  ur[n]  now 
contains  the  reflected  impulse  response  of  the  second  discontinuity.  We  can  ignore 
that  portion  of  the  reflection  response  caused  by  the  first  discontinuity,  and  record 
the  subsequent  data.  The  recording  window  is  defined  by  a  time  greater  than  the 
two-way  propagation  delay  time  to  the  first  discontinuity.  Next,  the  windowed  re¬ 
flection  response  is  deconvolved  across  the  first  discontinuity  transmission  impulse 
response  in  order  to  observe  the  second  discontinuity  reflected  impulse  response. 
Similarly,  another  AR  model  is  estimated  from  this  reflected  impulse  response,  and 
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then  we  repeat  the  procedure.  Practically,  the  excitation  signal  is  bandlimited  by 
the  instrumentation,  and  this  prevents  the  input  from  being  an  ideal  impulse. 


I  =  o 

vsM  j 


Inverse  Filter 


v0  tn] 

I 


V 

1 

T,(z) 

L 

1st  Discontinuity 


1  “4 


2nd  Discontinuity 


Figure  5.6  Layer-probing  of  a  unknown  discontinuity 

Also,  the  microstrip  dispersion  causes  pulse  spreading.  As  a  consequence,  when  two 
discontinuities  are  separated  by  less  than  the  time-dispersive  spread  for  the  propa¬ 
gating  pulse,  multiple  reflections  may  introduce  destructive  interference.  Therefore, 
when  the  reflected  impulse  response  is  significantly  distorted  by  destructive  inter¬ 
ference,  the  inverse  scattering  problem  cannot  be  solved. 

We  will  attempt  to  overcome  this  problem  by  measuring  the  output  trans¬ 
mission  response  resulting  from  the  last  impulsive  probing  signal.  The  remaining 
unknown  discontinuities  will  be  described  by  an  M-th  order  AR  model.  A  technique 
is  required  for  estimating  an  appropriate  model  order.  It  is  generally  observed  that 
M<N  where  N  is  the  data  segment  length.  The  best  order  is  typically  within  the 
range  from  N/20  to  N/5,  but  for  particularly  short  data  segments  the  best  order 
can  be  as  large  as  N/3  or  even  N/2  [Ref.  48].  Consequently,  we  do  have  a  priori 
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knowledge  about  the  order  by  assuming  the  model  order  and  the  detected  num¬ 
ber  of  discontinuities  are  equal.  As  discussed  earlier,  the  layer  reflection  coefficient 
algorithm  can  estimate  the  number  of  discontinuities  within  the  unknown  structure. 

The  transmission  response  now  reveals  sufficient  transient  information  to 
model  the  remaining  unknown  sections.  Using  the  WLS  algorithm,  the  transmission 
impulse  response  is  modeled  by  an  ARM  A  transfer  function  T(z).  Suppose  for 
simplicity,  and  not  too  much  loss  of  generality,  that  we  let  the  MA  model  (numerator 
term)  become  a  gain  constant,  so  that 

no  =  —  (5.65) 

and  the  transfer  function  can  be  rewritten  as 

M 

no  =  n  nw  (5.66) 

r=l 

where 


and  aT  is  the  AR  model  parameter.  As  can  be  seen,  the  individual  Tr(z)  functions 
are  lst-order  AR  transfer  functions  which  model  each  discontinuity.  Simulation 
results  are  performed  in  Chapter  VI  to  demo  '.strate  the  performance  of  the  layer¬ 
probing  algorithm. 

3.  Model  Attenuation  Factor  Estimation 

The  total  attenuation  factor  can  be  estimated  by  considering  the  energy  loss 
of  the  microstrip  structure.  An  incident  transient  pulse  delivers  a  finite  amount 
of  energy  to  the  microstrip  transmission  line.  Discontinuities  in  multi-section  mi¬ 
crostrip  structure  will  introduce  multiple  energy  reflections.  We  can  only  measure 
the  reflect  1  and  transmitted  voltage  waveforms  with  respect  to  the  input/output 


•r 
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measurement  ports.  Thus,  we  will  treat  the  lossy  multi-section  structure  as  a  sin¬ 
gle  cascaded  microstrip  step  discontinuity.  Referring  to  Figure  5.7,  the  cumulative 
energy  loss  is  given  by 

Eloss  =  e[,]  -  E{2t]  -  e[t)  (5.67) 

where 

E is  the  cumulative  reflected  energy  measured  at  the  input  port  (1), 

E^  is  the  cumulative  transmitted  energy  measured  at  the  output  port  (2), 

and 

e\'^  is  the  cumulative  incident  energy  measured  at  the  input  port  (1). 

The  cumulative  energy  values  are  calculated  by  [Ref.  31:  p.  34] 

N 

£  =  ^]u[n]|2  (5.68) 

»=0 


where  u[n]  is  the  sampled  voltage  measured  at  the  respective  input/output  pox-ts. 
Next,  the  energy  loss  is  described  in  terms  of  an  imaginary  discontinuity  boundary 
(D)  as 

Eloss  =  (4‘  -  4’)  +  (4*  -  4’)  +  (4*  -  4’)  (5-69) 

where 

r(d  _  r’(«),-2ai  /j 

^ D  ~  t 

r(r)  _  p(T)  2a\  l\ 
l  D  ~  c 

17(0  _  rlO  jai  In 

I,D  -  &  1  e 


Assuming  a  symmetrical  microstrip  geometry  with  respect  to  the  boundary,  we  will 
assume  Qj  =  ct>  and  /j  =  /2.  Equating  (5.67)  and  (5.69)  yields  the  half-symmetrical 
loss  (  in  nepers  )  as 


O]  h  = 


-L 
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E\r)  +  E\ 


(O' 


E\ 


(«) 


(Np) 


(5.70) 
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The  half-symmetrical  attenuation  factor  will  be  used  to  design  the  resistive  n- 
networks  in  Figure  4.5. 
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Figure  5.7  Cascaded  step  energy  flow  model 

This  chapter  presented  the  use  of  several  ARMA/AR  model  parameter  esti¬ 
mation  algorithms  for  the  characterization  and  transient  response  modeling  of  mi¬ 
crostrip  transmission  lines.  Both  deterministic  and  stochastic  approaches  were  dis¬ 
cussed.  Newly  developed  algorithms  to  improve  the  identification  of  lossy  microstrip 
lines  consisting  of  multiple  discontinuities  were  also  derived.  These  algorithms  were 
the  layer  reflection  coefficient  estimation  algorithm,  and  the  layer-probing  algorithm. 
The  purpose  of  these  algorithms  is  to  estimate  the  appropriate  ARMA/AR  model 
parameters  from  available  time-domain  measurements.  In  Chapter  VI,  the  algo¬ 
rithms’  performance  will  be  studied  using  simulated  and  experimentally  measured 
data. 
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VI.  EXPERIMENTAL  AND  SIMULATION  RESULTS 


The  purpose  of  this  chapter  is  to  verify  the  proposed  microstrip  models  using 
both  experimental  and  simulation  results.  The  lumped- distributed  equivalent  cir¬ 
cuit  and  parametric  models  presented  are  now  simulated  for  an  integrated  circuit 
(IC)  interconnection,  a  cascaded  microstrip  step  discontinuity,  and  a  multi-section 
microstrip  step  discontinuity.  The  microstrip  models  are  verified  by  comparing  pi¬ 
cosecond  time-domain  measurements  with  transient  analysis  and  parametric  model 
simulation  results.  Finally,  the  performance  of  the  ARMA/AR  parameter  estima¬ 
tion  algorithms  is  studied  using  simulated  and  measured  data. 

A.  INTEGRATED  CIRCUIT  INTERCONNECTION  RESULTS 

1.  Optoelectronic  Measurement  Techniques 

D.H  Auston  and  A.M.  Johnson  of  Bell  Laboratories  first  demonstrated  the 
generation  of  picosecond  pulses  on  silicon  substrates  in  1975  by  using  optically  trig¬ 
gered  gating  [Ref.  49,50].  They  used  quasimetallic  photoconductivity  produced  by 
the  absorption  of  short  optical  pulses  in  silicon  to  both  turn  on  and  off  a  switching 
gap.  Figure  6.1  shows  an  optoelectronic  switch  using  this  technique.  The  microstrip 
transmission  line  is  fabricated  on  high  resistivity  (p  —  104fl  —  cm)  silicon  and  con¬ 
tains  a  gap  which  forms  the  switch.  Optical  pulses  generating  mobile  carriers  within 
the  silicon  produce  the  switching  action.  The  work  on  silicon  has  been  extended  by 
W.R.  Eisenstadt  and  R.B.  Hammond  to  produce  picosecond  optoelectronic  switches 
on  bulk  silicon  [Ref.  51].  Their  effort  was  focused  towards  the  production  of  high 
speed  photoconductors  which  can  be  integrated  along  side  other  devices  on  an  in¬ 
tegrated  circuit  using  standard  VLSI  fabrication  techniques. 
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The  Photoconductor  Circuit  Element  (PCE)  is  the  best  laser  driven  optoelectronic 
transducer  for  high-speed  sampling  applications  to  integrate  on  a  silicon  substrate 
[Ref.  52].  The  primary  goal  of  the  PCE  is  to  achieve  an  ideal  impulsive  sampling 
gate. 


Figure  6.1  Diagram  of  D.H.  Auston's  picosecond  photoconductor 

2.  Picosecond  Pulse  Propagation  Measurements 

D.  Bowman  [Ref.  53:  pp.  82-94]  performed  a  series  of  picosecond  time- 
domain  measurements  using  a  polysilicon  PCE  test  structure.  The  structure  shown 
in  Figure  6.2  illustrates  the  IC  fabrication  of  the  PCE  sampler  on  the  microstrip 
transmission  line.  The  length  of  microstrip  beyond  the  last  PCE  sampler  prevents 
any  reflections  from  the  open  stub  end  from  returning  in  time  to  interfere  with 
the  measurement.  The  experimental  sampling  system  uses  a  colliding  pulse  mode- 
locked  (CPM)  ring  dye  laser  with  a  nominal  200  femtosecond  optical  pulse  to  excite 
a  3.5  GHz  bandwidth  eiectrical  pulse.  The  half  power  bandwidth  of  the  electrical 
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electrical  pulse  is  determined  from  the  magnitude  spectrum  using  the  Fast  Fourier 
transform.  Table  7  lists  the  physical  properties  of  the  microstrip  test  structure. 


Figure  6.2  Polysilicon  PCE  microstrip  test  structure 

a.  Sampled  Waveform  Measurements 

A  transient  pulse  was  measured  at  different  distances  along  a  100  //m 
wide  microstrip  at  a  sampling  interval  of  1.25  ps.  Figure  6.3  shows  the  measured 
waveforms  with  magnitudes  normalized  to  the  peak  of  the  first  sampler.  Successive 
curves  are  the  sampled  measurements  at  the  second  (500  /im),  third  (1000  /xm),  and 
fourth  (2000  /xm)  PCE  sampler. 
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NORMALIZED 


Table  7.  IC  MICROSTRIP  PROPERTIES 


Material  Properties 

Value 

Characteristic  Impedance  (Zq) 

73.24  Q 

Dielectric  Constant 

(M 

11.7 

Substrate  Thickness 

(h) 

350.5  fxm 

Conductor  Thickness 

(0 

1.0  fim 

Conductor  Width 

(w) 

100.0  /im 

Substrate  Resistivity  (Si) 

(Ps) 

14.5  fl-cm 

Conductor  Resistivity  (Al)  ( pc ) 

0.262e-05  fl-cm 

Conductor  Length 

(0 

500  /im 

Figure  6.3  High  frequency  waveform  propagation  measurements 
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The  measured  responses  show  the  expected  dispersive  propagation  with  high  losses. 
This  is  evident  from  the  increasing  rise  times  of  the  pulses  (pulse  spreading)  and 
smaller  peak  magnitudes  with  increasing  propagation  distance  [Ref.  53:  pp:  84-85]. 


b.  Propagation  Velocity 

The  rate  at  which  a  pulse  propagates  down  an  IC  interconnection  is  a  very 
important  parameter,  particularly  for  circuits  involving  critical  timing  applications 
or  pulse  synchronization.  At  the  high  frequencies  present  with  picosecond  pulses, 
propagation  velocity  varies  with  frequency  and  affects  the  shape  of  the  pulse. 

The  relative  time  at  which  the  peak  magnitude  of  the  pulse  occurs  can  be 
used  to  determine  the  velocity.  The  group  velocity  and  not  the  phase  velocity  of  the 
pulse  is  actually  measured  by  observing  the  relative  delays  of  the  peak  magnitude. 
Figure  6.4  reports  the  results  of  the  measured  propagation  delays  at  the  various 
time  differences  between  the  peak  magnitudes  for  each  of  the  three  samplers.  The 
theoretical  propagation  delay  is  calculated  by  [Ref  18:  p.  62] 


,  VQ71 

U  =  - 

c 


(6.1) 


for  the  different  lengths  of  the  microstrip.  What  is  interesting  about  the  results 
is  the  almost  identical  slope  for  both  the  measured  and  theoretical  delay.  This 
indicates  a  relatively  constant  velocity  of  w  112-114  ^m/ps. 

c.  Propagation  Loss 

The  decrease  in  the  peak  magnitude  of  the  propagating  pulse  with  dist  ance 
seen  in  Figure  6.3  is  evidence  of  the  significant  losses  present  on  silicon  substrates. 
Figure  6.5  shows  the  loss  in  decibels  as  a  function  of  propagation  distance.  The 
theoretical  loss  combines  the  numerical  results  for  conductor  loss  (2.8b)  and  dielec¬ 
tric  loss  (2.10)  at  three  microstrip  distances.  The  measured  propagation  losses  were 
obtained  from  relative  cumulative  energy  results.  The  theoretical  values  reveal  a 
loss  of  2.4  dB/mm  while  the  measurements  indicate  a  3.6  dB/mm  loss  is  actually 
occurring.  A  possible  explanation  for  the  discrepancy  between  the  theoretical  loss 
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and  the  measured  loss  is  that  (2.8b)  yields  a  somewhat  low  result  because  surface 
roughness  is  not  considered,  and  the  effective  microstrip  permittivity  value  used  in 
(2.10)  actually  varies  with  frequency  [Ref.  18:  p.90]. 
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Figure  6.4  Pulse  propagation  delay  on  IC  interconnection 

3.  Equivalent  Circuit  Simulations 

The  lumped-parameter  and  lumped-distributed  equivalent  circuits,  derived 
in  Chapter  II,  are  developed  for  a  500  /im  length  of  microstrip  in  the  PCE  test 
structure.  Numerical  results  for  the  circuit  model  axe  presented  in  Table  8.  Several 
comments  can  be  made  about  these  theoretical  calculations.  The  static  effective  mi¬ 
crostrip  permittivity  significantly  influences  all  subsequent  calculations.  It  assumes 
a  TEM  mode  of  propagation  on  the  microstrip  which  is  frequency  independent. 
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In  the  equivalent  circuit  model  the  lumped-elements  axe  evaluated  at  a  maximum 
frequency  of  3.5  GHz  which  corresponds  to  the  bandwidth  of  the  propagating  signal. 
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Figure  6.5  Pulse  propagation  loss  on  IC  interconnection 


At  this  frequency  the  physical  microstrip  length  should  be  less  than  4  mm  («  Aff/8). 
Consequently,  the  equivalent  circuit  models  may  lose  their  accuracy  if  significant 
higher  order  modes  are  excited  by  microstrip  discontinuities.  Figure  6.6  shows 
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Table  8.  EQUIVALENT  CIRCUIT  THEORETICAL  RESULTS 


Circuit  Parameters 

Equation 

Value 

Maximum  Frequency 

(/) 

3.5  GHz 

Static  Effective  Permittivity  ( 

£e//) 

(2.2) 

7.082 

Effective  Permittivity  (ce//(/)) 

(2.6) 

7.120 

Characteristic  Impedance 

(Zo) 

(2.3a) 

75.22  Q 

Guide  Wavelength 

(Af) 

(2.19) 

32.12  mm 

Shunt  Capacitance 

(Cr) 

(2.18a) 

0.03  pF 

Series  Inductance 

(Lr) 

(2.18b) 

0.33  nH 

Propagation  Loss 

Conductor  Loss 

(<*c) 

(2.8b) 

0.009  dB 

Dielectric  Loss 

M 

(2.10) 

1.204  dB 

Total  Loss 

1.213  dB 

the  schematic  diagrams  of  the  equivalent  circuit  models. 
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0.332  nH 


74.23H 


(a) 

Figure  6.6  (a)  Lumped-element  equivalent  circuit 

In  modeling  the  IC  interconnection  using  equivalent  circuits,  it  is  assumed  that 
impedance  matching  networks  are  used  to  prevent  reflections.  In  Figure  6.6(a)  a 


78 


propagation  loss  of  1.213  dB  is  modeled  by  a  resistive  7r-network  and  a  propagation 
delay  of  4.46  ps  is  described  by  the  reactive  LC  network.  Transient  analysis  simu- 


vs(t) 


74.23H 

vW — °- 


io.4n 

-VW- 


1.065  K12 


0.06  pF 


74.23Q 


(b) 

Figure  6.6(b)  Lumped-distributed  equivalent  circuit 

lations  of  the  circuit  models  are  performed  using  the  PSPICE  electrical  circuit  simu¬ 
lation  program  [Ref.  17].  The  PSPICE  network  listings  for  both  equivalent  circuits 
are  given  in  Appendix  B.  A  piecewise-linear  approximation  of  the  normalized  refer¬ 
ence  pulse,  measured  at  the  first  PCE  sampler,  acts  as  the  PSPICE  voltage  source. 
Waveform  comparisons  between  the  second  PCE  sampler  measurement  and  the 
circuit  model  simulations  are  shown  in  Figure  6.7.  The  general  shape  of  the  simu¬ 
lations  are  representative  of  the  measured  data;  however,  there  are  some  significant 
discrepancies.  Two  possible  causes  are  suggested.  First,  on  the  decaying  transient 
the  measurement  contains  a  second  peak  caused  by  a  reflected  signal  from  the  back¬ 
plane  of  the  substrate.  This  phenomena  is  not  modeled  by  the  circuit  model  because 
of  its  complex  wave  nature  [Ref.  53:  p.90j.  Second,  the  equivalent  T-network  circuit 
model  simulation  indicates  a  decreased  peak  magnitude.  This  observation  can  be 
explained  as  energy  lost  by  a  reflected  voltage  caused  by  the  shunt  capacitance. 
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Next,  we  derive  an  equivalent  circuit  model  for  the  1000/im  length  of  microstrip. 
Using  the  same  reactive  7r-network  of  Figure  6.6(a),  we  calculate  the  new  shunt 
capacitances  and  series  inductance  as  0.06  pF  and  0.663  nH,  respectively.  The 


Figure  6.7  Waveform  simulations  of  the  equivalent  circuit  models 

measured  propagation  loss  of  3.891  dB  is  used  to  evaluate  the  resistive  7r-network. 
Figure  6.8(a)  compares  the  PSPICE  transient  analysis  simulation  with  the  measured 
response.  The  simulated  results  tend  to  underestimate  the  rise  time  of  the  measured 
pulse. 
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This  error  may  be  attributed  to  the  breakdown  in  the  lumped- element  equivalent 
circuit  assumption  as  the  physical  length  of  the  microstrip  increases.  Therefore,  we 
can  improve  the  model  by  cascading  two  reactive  7r-network  sections.  A  PSPICE 
simulation  was  performed  using  the  higher  order  circuit  model. 


Figure  6.8(a)  Single  LC  network  simulation  at  the  1000/xm  PCE  sampler 

Figure  6.8(b)  shows  a  significant  improvement  in  the  waveform  agreement,  especially 
in  the  rise  times.  However,  as  the  number  of  reactive  lumped-elements  increase,  the 
model  order  also  increases. 
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NORMALIZED  AMPLITUDE  (V) 


Then  the  mathematical  relationships  among  the  model  parameters  and  the  mi¬ 
crostrip  physical  properties  will  become  more  complex. 


Figure  6.8(b)  Cascaded  LC  network  simulation  at  the  1000/im  PCE  sampler 
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4.  Parametric  Model  Simulations 


The  following  ARMA  digital  transfer  function  is  computed  for  the  500  /zm 
microstrip  test  section  using  (3.10)  and  (3.11)  at  a  sampling  interval  of  1.25  ps: 

0.1172 


y-i 


Harma{z)  = 


1  -  1.4513  2"1  +  0.57  2 


—2  ' 


(6.2) 


The  normalized  reference  pulse,  measured  at  the  first  PCE  sampler,  is  attenuated 
by  1.213  dB  and  then  iterated  through  a  difference  equation  of  (6.2).  Similarly,  an 
AR  digital  transfer  function  is  evaluated  using  (3.19)  as 


HAr(z) 


0.4295 

1  -  0.5705  z-1' 


(6.3) 


The  normalized  reference  pulse  is  again  attenuated  by  1.213  dB  and  iterated  through 
a  difference  equation  of  (6.3).  Figure  6.9  shows  the  ARMA  and  AR  model  simulation 
results  as  compared  with  the  measured  pulse. 


Figure  6.9(a)  lst-order  ARMA  parametric  model  simulation 
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Simulated  pulses  are  shown  by  solid  lines  while  the  measurements  are  indicated 
with  dash  lines.  Both  simulated  waveforms  have  properly  accounted  for  the  loss 
mechanisms  on  the  500/xm  microstrip  section;  however,  the  2nd-order  ARMA  model 
result  shews  an  outstanding  waveform  agreement  as  compared  with  the  Ist-order 
AR  model.  In  the  next  section  we  will  focus  on  the  model  parameter  estimation 
algorithms  using  both  simulated  and  measured  data. 


Figure  6.9(b)  lst-order  AR  parametric  model  simulation 


5.  Model  Parameter  Estimation  Results 
a.  WLS  Algorithm 

The  weighted  least  squares  (WLS)  algorithm  of  Chapter  V  can  be  used  to 
estimate  model  parameters  from  a  given  sample  impulse  response.  First,  we  must 
obtain  an  accurate  impulse  response  from  the  available  input/output  measurements. 
The  impulse  response  was  obtained  by  deconvolving  the  output  data,  measured  at 
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the  second  PCE  sampler,  with  the  measured  normalized  input  data.  The  original 
measured  data,  sampled  at  1.25  ps,  was  decimated  to  a  sampling  interval  of  5 
ps  to  improve  the  frequency  resolution  of  the  FFT  during  deconvolution.  Next, 
the  resulting  sample  impulse  response  was  applied  to  the  WLS  algorithm.  The 
estimated  ARMA  model  parameters  are  listed  in  Table  9.  Figure  6.10  compares  the 

Table  9.  ARMA  PARAMETERS  USING  WLS  ALGORITHM 

Theoretical 

0.657 
0.285 
-0.106 


Estimated 

0.688 

0.187 

-0.107 


Parameter 

bx 

ax 

Q2 


estimated  impulse  response  to  the  measured  and  theoretical  impulse  responses  at  a 
5  ps  sampling  interval. 


Figure  6.10  Comparison  of  ARMA  model  sample  impulse  responses 
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b.  Stochastic  Methods 


ARMA  and  AR  sample  impulse  responses  were  simulated  using  a  PSPICE 
transient  analysis  of  the  equivalent  circuit  models  at  a  sampling  interval  of  1.25  ps. 
Stochastic  output  processes  are  generated  using  the  convolution  summation 

N 

yW  =  £  i>[A']/i[n  —  k]  (6.4) 

k= o 


where  r>[n]  is  an  input  Gaussian  white  noise  sequence  and  h[ri]  is  the  simulated 
sample  impulse  response.  The  white  noise  sequence  was  obtained  using  a  pseudo¬ 
random  number  generator  program.  The  output  process  t/[n],  having  a  record  length 
of  1800  samples,  and  the  original  white  noise  sequence  were  applied  to  the  M-R 
algorithm.  The  ARMA  digital  filter  coefficients  were  estimated  for  a  (2,1)  order 
model.  Table  10  compares  these  estimates  to  the  theoretical  parameter  values  at  a 
1.25  ps  sampling  interval. 

Table  10.  ARMA  PARAMETERS  USING  M-R  ALGORITHM 


Parameter 

Theoretical 

Estimated 

bi 

0.1172 

0.1137 

1.4513 

1.4796 

°2 

-0.5705 

-0.5905 

The  estimated  results  compare  favorably  with  the  theoretical  model  parameter  val¬ 
ues.  There  is,  however,  a  limitation  in  using  the  stochastic  approach.  First,  the 
M-R  algorithm  requires  input/output  data  to  be  simultaneously  sampled  to  provide 
accurate  ARMA  parameter  estimates.  Second,  microwave  measurement  techniques 
are  usually  based  on  frequency- domain  scattering  parameters,  and  high-speed  pi¬ 
cosecond  sampling  techniques  are  not  commercially  available. 

The  Schur  algorithm  is  used  to  estimate  a  lst-order  AR  model.  Initially,  a 
lst-order  AR  sample  impulse  response  was  simulated  using  PSPICE  at  a  sampling 
interval  of  1.25  ps.  Using  equation  (5.45),  a  sample  autocorrelation  sequence  is 
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computed  from  the  impulse  response  data.  Then  the  Schur  algorithm  is  used  to 
solve  the  AR  reflection  coefficients,  and  the  normalized  AR  transfer  function  yields 

H •  <6-5> 

Since  the  impulse  response  was  simulated  using  a  lst-order  AR  model  a  single  Schur 
reflection  coefficient  should  be  expected. 

c.  Elmore  Method 

Network  functions  are  now  estimated  from  rise  and  delay  time  measure¬ 
ments.  Figure  6.11  describes  two  consecutive  PCE  sampler  voltage  measurements 
denoted  as  (t)  and  (t).  A  two-port  network  describes  the  500/im  length  of 
microstrip. 


Figure  6.11  Cascaded  network  functions 

The  Elmore  rise  time  ( Tr )  is  defined  as  the  reciprocal  of  the  slope  of  the  tangent 
drawn  to  the  response  curve  at  its  half-magnitude  point  [Ref.  46].  The  delay  times 
are  estimated  at  consecutive  peak  or  half-magnitude  points.  Referring  to  Figure 
6.12,  a  delay  time  of  4.41  ps  is  obtained  by  averaging  the  peak  magnitude  delay 
(4.48  ps)  and  the  half-magnitude  (4.33  ps)  delay  times.  A  rise  time  of  2.06  ps  is 
computed  using  (5.34)  for  Tr 2  =  7.73  ps  and  TV3  =  8.00  ps. 
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A  2nd-order  network  function  is  evaluated  as 


H(s) 


99.3836  x  1021 

s2  +  438.282  x  109s  +  99.3836  x  1021 ' 


(6.6) 


The  characteristic  impedance  of  the  microstrip  must  be  determined  prior  to  solv¬ 
ing  for  uiq  and  £  in  (3.5).  The  estimation  of  the  characteristic  impedance  will  be 
addressed  in  the  next  section.  Finally,  the  ARMA  model  parameters  are  calculated 
at  the  sampling  interval  of  5  ps  using  equations  (3.10)  and  (3.11). 


Figure  6.12  Elmore  delay  and  rise  time  graphical  estimates 

Equation  (5.35)  solves  the  lst-order  AR  filter  coefficient  in  terms  of  the  delay  and 
rise  times.  Both  the  estimated  and  theoretical  ARMA/AR  model  parameters  are 
compared  in  Table  11  at  a  1.25  ps  sampling  interval. 
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The  Elmore  estimates  for  the  AR  model  parameters,  including  the  denominator  AR 
filter  coefficients  of  the  ARMA  model,  all  show  a  slight  increase  over  their  theoretical 


Table  11.  MODEL  PARAMETERS  USING  ELMORE 


ARMA  Parameter 

Theoretical 

Estimated 

bi 

0.1172 

0.1165 

°l 

1.4513 

1.4602 

°2 

-0.5705 

-0.5782 

AR  Parameter 

a 

0.5705 

0.5782 

values.  This  result  will  explain  the  increase  in  effective  microstrip  permittivity 
estimates  presented  in  the  next  section. 

d.  Estimation  of  Physical  Microstrip  Properties 

The  characteristic  impedance  at  each  PCE  sampler  is  computed  using 
(2.29)  and  (2.30).  Table  12  lists  the  calculated  characteristic  impedance  estimated 
at  each  PCE  sampler.  The  average  impedance  of  74.32  fl  compares  very  favorably 
to  the  actual  microstrip  characteristic  impedance  of  74.23  Q. 

Table  12.  CHARACTERSITIC  IMPEDANCE  ESTIMATES 


PCE  Sampler  No. 

Zo 

(1) 

72.92  U 

(2) 

73.03  n 

(3) 

77.06  n 

Average  Value 

74.34  ft 

The  effective  microstrip  permittivity  (ee//)  is  directly  calculated  from  the 
model  parameter  estimates  using  (3.22)  and  (3.25).  Several  theoretical  expressions 
have  been  cited  in  the  literature  for  the  effective  microstrip  permittivity.  Here,  we 
have  computed  a  value  of  7.082  using  equation  (2.2).  However,  other  closed-form 
expressions  will  produce  different  results.  For  example,  another  expression  for  the 
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effective  microstrip  permittivity,  given  the  characteristic  impedance,  is  [Ref.  IS: 

P-44] 


e<//  = 


£r  +  1 


1  + 


29.98 


Zq  V  Cr  +  1 


f-r  ~  1 
tr  +  1 


■  i  1  .  4 
In  -  +  —In  - 

2  fr  7T 


(6.7) 


£e// 


=  7.186 


where  Zq  =  74.23  Q  and  tT  =  11.7  for  silicon  substrate.  Table  13  lists  the  effective 
microstrip  permittivity  values  that  were  computed  using  the  ARMA/AR  model 
parameters  obtained  from  the  different  estimation  algorithms. 


Table  13.  EFFECTIVE  PERMITTIVITY  ESTIMATES 


Estimator 

ARMA 

AR 

WLS 

7.210 

Schur 

7.087 

Elmore 

7.511 

7.496 

These  results  indicate  that  the  best  model  parameter  estimators  are  the  WLS  and 
Schur  algorithms.  These  algorithms  provide  optimal  estimates  in  a  least  squares 
sense,  while  the  Elmore  method  was  extremely  dependent  upon  obtaining  accurate 
rise  and  delay  time  measurements.  Precautions  should  be  made  whenever  we  com¬ 
pare  the  estimated  results  to  the  theoretical  values.  For  example,  the  accuracy  of 
the  theoretical  expressions  is  dependent  upon  the  shape  ratio  w/h  range.  In  all 
cases  the  shape  ratio  will  be  accurate  to  dhl  percent.  For  narrow  lines  (w/h  <  1.3), 
the  effective  microstrip  permittivity  has  an  error  range  +0.5-0.0  percent.  When 
calculated  using  (6.7),  ceff  is  accurate  to  ±1  percent  [Ref.  18:  pp.  45-46].  Ad¬ 
ditionally,  the  effective  microstrip  permittivity  is  in  fact  frequency  dependent,  and 
the  values  tend  to  be  slightly  higher  than  those  given  by  (2.2)  and  (6.7). 
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B.  CASCADED  MICROSTRIP  STEP  DISCONTINUITY  RESULTS 

The  equivalent  circuit  and  parametric  models  of  Chapter  IV  are  validated  by 
comparing  model  simulations  with  experimental  measurements.  Figure  6.13  de¬ 
scribes  a  cascaded  microstrip  step  discontinuity  structure  that  was  fabricated  on 
G-10  epoxy  dielectric  material.  The  physical  dimensions  of  the  cascaded  step  dis¬ 
continuity  were  restricted  by  the  available  photolithographic  equipment  and  the 
dielectric  substrate.  A  half  power  bandwidth  of  400  MHz  for  the  measured  refer¬ 
ence  pulse  was  selected  as  the  maximum  design  frequency.  This  frequency  is  used  in 
the  theoretical  calculations.  Table  14  summarizes  the  theoretical  microstrip  results 
obtained  from  the  closed-form  expressions. 


Figure  6.13  Cascaded  microstrip  step  discontinuity 

1.  Validation  of  Equivalent  Circuit  Models 

The  lumped-element  values  in  Table  14  are  associated  with  the  equivalent 
circuit  models  of  Figures  4.3(c)  and  4.4(b).  A  transient  analysis  simulation,  using 
PSPICE,  is  performed  for  each  equivalent  circuit.  A  piecewise-linear  approximation 
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of  the  measured  input  pulse  was  used  in  the  transient  analysis  simulations.  A 
half-symmetric  attenuation  factor  of  0.0348  nepers  is  computed  from  the  measured 
incident,  reflected,  and  transmitted  cumulative  energy.  Separate  resistive  at  tenuator 
7r-networks  arc  designed  to  model  each  half-symmetric  loss.  The  shunt  resistors 
(2.316  kfl)  are  given  by  [Ref.  37:  pp. 186-192]: 

«>  =  z«  (£rr)  <c-8»> 

and  the  series  resistor  (1.41  fl)  by 

l> 2  =  Z.  (Kjr)  (6.86) 

where  the  loss  in  decibels  is  201og10  N.  Transmission  and  reflection  output  re¬ 
sponses  are  simulated.  These  simulations  arc  compared  to  the  measured  responses 
in  Figures  6.14(a)  and  6.14(b),  respectively.  The  measured  transmission 
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Figure  6.14(a)  PSPICE  simulated  transmission  response 
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and  reflection  responses  were  obtained  using  the  inverse  FFT  of  the  measured  scat¬ 
tering  "..rameters. 

2.  Parametric  Model  Simulations 

Using  equations  (4.14),  and  (4.16)  —  (4.18),  yields  a  3rd-order  ARMA  trans¬ 
fer  function 

^  x  6.776  x  109  -  7.295  x  lO9^1  +  4.323  x  109c~2 

73{r)~  i  _i.702c-i  +  1.59643-2  -  0.762c-3  (6'9) 

at  a  sampling  interval  of  40  ps.  A  sample  impulse  response  is  generated  using  the 
difference  equation  of  (6.9).  The  sampled  incident  pulse  is  then  convolved  with  the 
sample  impulse  response  data  to  produce  an  output  transmission  response. 


Table  14.  CASCADED  STEP  THEORETICAL  RESULTS 


Microstrip  Parameter 


Equation 


Wide  Width  (3  cm)  Segment 

Static  Effective  Permittivity 

(Ceff) 

Effective  Permittivity 

(Uff(f)) 

Characteristic  Impedance 

(Zn) 

Guide  Wavelength 

(V>i) 

Shunt  Capacitance 

(Cs) 

End  Inductance 

(Lc) 

Value 


(2.2) 


(2.6a) 


(2.3a) 


(2.19) 


(4.7) 


(4.8) 


3.953 


4.024 


8.44  D 


373.90  mm 


13.10  pF 


0.931  nH 


Narrow  Width  (0.437  cm)  Segment 


Static  Effective  Permittivity 

(2.2) 

3.354 

Effective  Permittivity 

(««//(/)) 

(2.6a) 

3.360 

Characteristic  Impedance 

(Z02) 

(2.3b) 

40.39  n 

Guide  Wavelength 

(^g02  ) 

(2.19) 

409.16  mm 

Series  Inductance 

(L,) 

(4.5) 

8.640  nH 

Step  Discontinuity  Circuit  Parameters 


Shunt  Capacitance  ( Cd ) 


Series  Inductance  ( L\ ) 


Series  Inductance  (Z-?) 


(Cd) 

(4.3a) 

0.559  pF 

(Li) 

(4.1a) 

0.054  nH 

{L2) 

(4.1b) 

0.235  nH 

The  normalized  transmission  response  is  compared  to  the  measured  output  in  Figure 
6.15.  Next,  the  lst-order  ARMA  reflection  and  transmission  transfer  functions  are 


Figure  G.  14(b)  PSPICE  simulated  reflection  response 


obtained  using  (4.20b)  and  (4.20c),  respectively,  as 

=  (■-;-!) 
'  1  —  0.492z-1 

and 

0^54(1+^ 

v  '  1  -  0.492z_1 


(G.lOu) 


(6.105) 


at  a  sampling  interval  of  200  ps. 


Figures  6.16(a)  and  6.16(b)  compare  the  ARMA  reflection  and  transmission  re¬ 
sponse  simulations  with  the  measured  responses.  The  results  show  both  the  lossless 


-0  JO  I  I  I  I  I  I  i—i-  I  I  I  -I  I  I  I  I  i  i  i  i  I  i  i  i  i  I  r-r  i  i  I  I  i  i  i  i  i  i  i  i  1  i  i  i  i  i  i  i" i~ >  | 
OjOE+OOO  1.0E-009  2.0C— 009  3.06-009  4.06-009  6j0E-O09 

TIME  (  SECOND  ) 


Figure  6.15  3rd-order  ARMA  model  transmission  response 

and  lossy  (0.0348  nepers)  simulated  output  waveforms.  When  the  attenuation  loss 
is  included,  the  waveform  agreement  improves,  especially  with  respect  to  the  peak 
amplitudes. 

Finally,  a  lst-order  AR  model  is  given  by  (4.21)  as 

Tar(z)  = - °'68  T  (6.11) 

ARK  '  1  —  0.507z-1  V 

at  a  sampling  interval  of  200  ps.  The  normalized  AR  transmission  response  is 
compared  to  the  measured  response  in  Figure  6.17. 
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Figure  6.16  lst-order  ARM  A  (a)  reflection  response  and  (b)  transmission  response 
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3.  Estimation  Results 

a.  Estimation  of  Reflection  Impulse  Response 

An  appropriate  parametric  model  for  the  cascaded  microstrip  step  dis¬ 
continuity  is  estimated  from  only  the  measured  incident  and  reflected  responses. 
An  estimate  of  the  reflection  sample  impulse  response  is  obtained  by  deconvolving 
the  measured  reflection  response  with  the  incident  pulse.  Figure  C.18  shows  the 
estimated  reflection  sample  impulse  response  at  a  sampling  interval  of  200  ps. 
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Figure  G.17  lst-order  AR  transmission  response 
b.  Estimation  of  Model  Parameters 

The  WLS  algorithm  estimates  the  lst-order  AR  model  parameter  of  equa¬ 
tion  (4.21)  using  the  previously  estimated  reflection  impulse  response.  The  time 
constant  ( C\Zq2 )  is  solved  as  291.51  ps  using  an  estimated  lst-order  AR  filter  co¬ 
efficient  of  0.5035. 


97 


Next  the  lst-order  ARMA  model  parameters  axe  calculated  by  substituting  this 
time  constant  value  into  equation  (4.20c).  Therefore,  the  unknown  transmission 
response  can  also  be  approximated  from  observing  only  the  incident  and  reflected 
data.  Finally,  the  ARMA  and  AR  model  parameter  estimates  are  given  in 


Figure  6.18  Estimated  reflection  sample  impulse  response 

Table  15.  This  approach  will  be  exploited  by  the  layer-probing  algorithm  in  the 
next  section  for  the  system  identification  of  a  multi-section  microstrip  structure 
[Ref.  54]. 
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Table  15.  MODEL  PARAMETER  ESTIMATES 


ARMA  Parameter 

Theoretical 

Estimated 

bo 

0.2538 

0.2554 

<21 

0.4924 

0.4892 

AR  Parameter 

a 

0.5065 

0.5035 

C.  MULTI-SECTION  STEP  DISCONTINUITY  RESULTS 

As  we  saw  earlier,  a  cascaded  step  microstrip  discontinuity  can  be  modeled  by 
an  equivalent  circuit  or  by  an  AR  parametric  model.  We  will  now  extend  these 
models  for  a  multi-section  microstrip  step  discontinuity. 

1.  Equivalent  Lumped-Distributed  Circuit  Model 

a.  Theoretical  Results 

The  capacitive  shunt  equivalent  circuit  model,  as  shown  in  Figure  4.4, 
will  be  used  to  characterize  each  cascaded  step  discontinuity.  The  high  impedance 
(narrow  microstrip  width)  lines  which  connect  each  discontinuity  are  modeled  as 
lossless  distributed  transmission-line  sections  using  (2.24).  The  propagation  delay 
time  associated  with  each  transmission  line  is  determined  from  the  physical  length 
and  the  frequency-dependent  effective  microstrip  permittivity.  Figure  6.19  describes 
the  multi-section  microstrip  structure  which  will  be  modeled  and  simulated  using 
PSPICE.  Three  different  line  width  dimensions  are  referred  to  by  numbered  line 
segments.  Specific  equivalent  circuit  numerical  results  are  fisted  in  Table  16  for 
each  fine  segment.  The  theoretical  calculations  use  a  maximum  design  frequency 
of  4  GHz  and  a  substrate  relative  dielectric  constant  er  of  4.3.  Frequency  disper¬ 
sion  is  incorporated  into  the  circuit  model  by  using  frequency  dependent  effective 
microstrip  permittivity  ce//(/)  and  characteristic  impedance  Zo(f )  values. 
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(6.12) 


The  characteristic  impedance  is  calculated  using  [Ref.  18:  pp.  81-82] 

376.6  h 

°U)  »ef;if)y/wU) 

where 

XI)  r  t  f  ~  tO 

weff(f)  =  w+ - 7TT2  ’ 

1+(i) 

f  = _ £ _ 

P  2 

and  the  effective  width  is 

376.7  h 
Weff  ~  Zo^JJ 


Figure  6.19  Multi-section  microstrip  discontinuity  test  structure 


The  theoretical  propagation  losses  axe  given  in  Table  17  for  each  line  segment.  The 
total  segment  loss  in  dB/Xg  is  the  sum  of  each  segment  dielectric  aj  and  conductor 
ac  loss,  given  by 

27.3  er  (ce/y  -  l)tan<5 

Qd  Ce//(Cr-!) 

0.072 
wZo(f) 

where  tan  8  is  the  dielectric  loss  tangent  and  /  is  in  GHz. 


(2.10) 

(2.86) 
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Table  16.  MULTI-SECTION  MICROSTRIP  THEORETICAL  RESULTS 


Microstrip  Parameter 

Equation 

Value 

Narrow  (0.3  cm)  Segment  No.  (1) 

Static  Effective  Permittivity 

V, 

(2.2) 

3.244 

Frequency-Dependent  Eff.  Permittivity 

<%(/) 

(2.6) 

3.358 

Static  Characteristic  Impedance 

z'1’ 

(2.3a) 

51.3  ci 

Frequency- Dependent  Impedance 

4°(/) 

(6.12) 

53.0  Cl 

Guide  Wavelength 

A<» 

(2.19) 

41.00  mm 

Wide  (3  cm)  Segment  No.  (2) 

Static  Effective  Permittivity 

(2.2) 

3.953 

Frequency- Dependent  Eff.  Permittivity 

#/(/) 

(2.6) 

4.100 

Static  Characteristic  Impedance 

Z«2) 

(2.3a) 

8.5  n 

Frequency-Dependent  Impedance 

42,(/) 

(6.12) 

9.5  Cl 

Guide  Wavelength 

x(2) 

(2.19) 

37.04  mm 

Shunt  Capacitance 

r(2) 

U5 

(4.7) 

2.14  PF 

Wide  (1.5  cm)  Segment  No.  (3) 

Static  Effective  Permittivity 

(2.2) 

3.755 

Frequency-Dependent  Eff.  Permittivity 

43A</> 

(2.6) 

3.922 

Static  Characteristic  Impedance 

43) 

(2.3a) 

15.7  Cl 

FVequency- Dependent  Impedance 

43\f> 

(6.12) 

17.5  Cl 

Guide  Wavelength 

\(3) 

Ag 

(2.19) 

37.87  mm 

Shunt  Capacitance 

fl3) 

(4.7) 

3.13  pF 

Step  Discontinuity  Circuit  Parameters 

Shunt  Capacitance 

Cf 

(4.3a) 

0.723  pF 

Shunt  Capacitance 

cf 

(4.3a) 

0.213  pF 

The  total  propagation  loss  (in  dB)  for  each  microstrip  line  is  modeled  in  the  circuit 
model  by  a  resistive  attenuator  x-networks.  Figure  6.20  shows  the  equivalent  circuit 

Table  17.  MULTI-SECTION  MICROSTRIP  THEORETICAL  LOSSES 


Microstrip  Attenuation 

Equation 

Value 

Narrow  (0.3  cm)  Segment  No.  (1) 

Dielectric  Loss 

(2.10) 

0.450  dB/A9 

Conductor  Loss 

c<‘» 

(2.8b) 

0.037  dB/Aff 

Total  Segment  Loss 

a<J> 

0.487  dB/A g 

Wide  (3  cm)  Segment  No.  (2) 

Dielectric  Loss 

(2) 

Qd 

(2.10) 

0.484  dB/A? 

Conductor  Loss 

EH 

(2.8b) 

0.019  dB/Aff 

Total  Segment  Loss 

(2) 
a\  ' 

0.503  dB/Aff 

Wide  (1.5  cm)  Segment  No.  (3) 

Dielectric  Loss 

(3) 

ad 

(2.10) 

0.477  dB/A g 

Conductor  Loss 

(2.8b) 

0.021  dB/A,, 

Total  Segment  Loss 

a<3> 

0.498  dB/A9 

model  for  the  microstrip  filter  of  Figure  6.19. 
b.  Simulation  Results 

The  microstrip  filter  was  fabricated  on  G-10  epoxy  dielectric  material. 
Scattering  parameter  measurements  were  taken  using  the  network  analyzer.  A 
piecewise-linear  approximation  of  the  measured  excitation  signal  is  used  in  the 
PSPICE  circuit  simulation.  The  transient  analysis  of  the  equivalent  circuit  model 
provides  both  reflection  and  transmission  responses.  Shown  in  Figure  6.21  is  a  com¬ 
parison  between  the  measured  and  simulated  reflection  responses.  The  amplitude  of 
the  simulated  response  indicates  that  the  theoretical  losses  are  slightly  greater  than 
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than  the  measured  losses.  The  magnitude  spectrum  of  the  reflection  responses  is 
shown  in  Figure  6.22. 


Figure  6.20  Multi-section  equivalent  circuit  model 

The  multi-section  microstrip  filter  is  also  simulated  by  PUFF,  a  computer- 
aided  design  (CAD)  program  for  microwave  integrated  circuits  [Ref.  55].  The  PUFF 
simulation  solves  the  scattering  parameters  for  the  microstrip  filter  described  in  Fig¬ 
ure  6.19.  The  inverse  FFT  is  performed  on  the  scattering  parameters  to  obtain  the 
reflected  and  transmitted  impulse  responses.  The  PUFF  program  input  excitation 
is  a  gaussian-shaped  unit  amplitude  pulse.  In  order  to  accurately  compare  both  the 
simulation  results,  a  piecewise- linear  approximation  of  the  PUFF  excitation  was 
used  in  the  PSPICE  transient  analysis.  Figures  6.23  and  6.24  compare  the  reflec¬ 
tion  and  transmission  responses,  respectively.  The  PUFF  simulation  model  does 
not  account  for  propagation  losses. 
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Figure  6.24  illustrates  the  effect  of  losses  in  the  equivalent  circuit  model.  The 
influt  .^e  of  loss  and  dispcision  on  the  multiple  reflected  signals  can  be  observed  in 
the  Figure  G.23  by  the  apparent  phase  shift  in  the  PSPICE  simulation. 


Figure  6.21  Comparison  of  PSPICE  reflection  response  waveforms 
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Figure  6.22  PSPICE  reflection  response  magnitude  spectrum 
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Figure  0.23  Comparison  of  simulated  reflection  responses 
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2.  Layer-Probing  Simulation  Results 

A  lumped-distributed  equivalent  circuit  is  developed  for  a  three-sectioned 
microstrip  step  discontinuity  structure.  A  PSPICE  transient  analysis  is  performed 


Figure  6.24  Simulated  transmission  responses 

using  the  equivalent  circuit  model.  Reflection  and  transmission  response  data  will  be 
applied  to  the  layer  reflection  coefficient  estimation  and  layer-probing  algorithms. 
The  initial  reflection  impulse  response  data  is  applied  to  a  Fortran  program  im¬ 
plementation  of  the  layer  reflection  coefficient  estimation  algorithm.  As  shown  in 
Figure  6.25,  a  sign  change  in  the  layer  reflection  coefficient  sequence  corresponds 
to  an  abrupt  width  change  on  the  microstrip  line.  An  estimated  group  velocity  of 
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2.53  mm/sample  is  computed  from  a  measured  propagation  time  delay  of  707.9  ps 
and  a  structure  length  of  112  mm.  Figure  6.26  shows  the  impedance  profile 


Figure  6.25  Multi-section  layer  reflection  coefficient  sequence 

computed  from  the  estimated  layer  reflection  coefficients.  The  accuracy  of  the 
impedance  profile  degrades  rapidly  beyond  the  first  discontinuty  because  of  loss  and 
dispersion.  Nevertheless,  we  axe  able  to  detect  the  relative  changes  in  the  impedance 
profile  and  then  estimate  the  number  of  discontinuities  within  the  structure.  This 
result  is  more  qualitative  than  quantitative.  The  detection  resolution  performance 
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of  the  algorithm  is  dependent  upon  the  sampling  interval  used.  Accurate  layer 
reflection  coefficients  are  obtained  only  when  the  sampling  interval  is  sufficiently 
small. 


Figure  6.26  Multi-section  impedance  profile 

The  WLS  algorithm  was  used  to  efficiently  compute  a  1st -order  AR  model 
from  the  reflected  impulse  response  shown  in  Figure  6.27.  The  layer-probing  algo¬ 
rithm  next  formulates  an  appropriate  inverse  filter  using  the  estimated  AR  model 
parameter.  A  new  excitation  signal  is  simulated  from  the  difference  equation  de¬ 
scription  of  the  first  inverse  filter  transfer  function.  We  create  a  piecewise- linear 
approximation  of  the  customized  input  and  initiate  a  second  transient  analysis  in 
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PSPICE  to  simulate  the  pulse  propagating  in  the  microstrip.  The  resulting  second 
reflection  response  is  windowed  to  eliminate  the  reflected  energy  caused  by  the  first 


Figure  6.27  Initial  reflection  response  waveform 

discontinuity.  The  windowed  reflection  response  is  shown  in  Figure  6.28.  The  reflec¬ 
tion  response  is  now  deconvolved  with  the  estimated  transmission  impulse  response 
of  the  first  discontinuity.  Figure  6.29  shows  the  result  of  the  deconvolution.  The 
WLS  algorithm  is  used  to  estimate  another  lst-order  AR  model  that  will  describe 
the  second  discontinuity.  The  previous  AR  model  parameter  estimates  axe  used  to 
design  a  2nd-order  inverse  filter.  Next,  we  excite  the  difference  equation  description 
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of  the  inverse  filter  using  the  original  input  pulse.  A  third  excitation  sequence  is  gen¬ 
erated  to  probe  the  microstrip.  Once  again  a  piecewise-linear  approximation  of  the 
third  customized  input  is  entered  in  PSPICE.  Figure  6.30  shows  the  transmission 


Figure  6.28  2nd  windowed  reflection  response 

response  simulated  by  the  multi-section  equivalent  circuit  model.  The  simulated 
transmission  response  immediately  reveals  the  transmitted  impulse  response  of  the 
third  discontinuity.  The  transmitted  impulse  response  is  identified  as  the  exponen¬ 
tially  decaying  transient  shown  in  Figure  6.30.  Finally,  a  third  AR  model  is 
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obtained  from  the  transmitted  impulse  response  using  the  WLS  algorithm.  The 
purpose  of  the  layer-probing  algorithm  was  to  resolve  the  reflected  and  transmitted 
impulse  response  for  each  cascaded  step  discontinuity.  This  was  accomplished 


Figure  6.29  2nd  reflected  impulse  response  after  deconvolution 

by  iterately  “probing'’  the  multi-section  structure  with  customized  signals  synthet¬ 
ically  produced  by  inverse  filtering.  Table  18  compares  the  estimated  AR  model 
parameters  with  the  theoretical  values  used  in  the  simulation. 


Ill 


3RD  LAYER  TRANSMISSION  RESPONSE  (V) 


The  multi-section  microstrip  structure  of  Figure  6.19  was  excited  by 
a  18.5  GHz  incident  pulse.  The  performance  of  the  layer  reflection  coefficient 

Table  18.  MULTI-SECTION  MODEL  PARAMETER  ESTIMATES 


Discontinuity  Parameter 

Theoretical 

Estimated 

ai  (T  =  2.0  ps) 

0.9829 

0.9816 

a2  (T  =  2.0  ps) 

0.9560 

0.9567 

fl3  (T  =  2.5  ps) 

0.9766 

0.9769 

Figure  6.30  Third  discontinuity  transmitted  impulse  response 


estimation  algorithm  is  studied  using  the  measured  reflection  response  shown  in 
Figure  6.31.  The  layer  reflection  coefficients  are  shown  in  Figure  6.32. 


Figure  6.31  Measured  reflection  response  due  to  the  18.5  GHz  pulse 

However,  the  estimated  reflection  coefficient  values  rapidly  decreased  in  amplitude 
after  the  30th  sample  value.  In  Figure  6.31,  we  can  observe  a  significant  initial 
reflection  which  is  followed  by  several  lower  amplitude  transients.  The  primary 
reflection  is  due  to  the  first  discontinuity,  while  the  subsequent  transients  are  dis¬ 
torted  by  the  multiple  reflections.  Next,  the  impedance  profile  is  computed  from 
the  estimated  reflection  coefficients  using  an  initial  impedance  value  of  50Q. 
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Referring  to  Figure  6.33,  a  characteristic  impedance  of  120  is  estimated  for  the  first 
step  discontinuity  as  compared  with  a  theoretical  value  of  9.5 0  at  4  GHz. 


Figure  6.32  Layer  reflection  coefficient  estimate  using  measured  data 

The  impedance  profile  estimates  a  500  microstrip  line  following  the  initial  disconti¬ 
nuity.  This  narrow  impedance  estimate  compares  favorably  to  the  theoretical  value 
of  530. 

The  simulation  and  measured  results  presented  in  this  section  have  demon¬ 
strated  the  performance  of  the  layer-probing  algorithm.  However,  the  algorithm 
assumes  that  the  multiple  discontinuities  can  be  accurately  modeled  by  cascading 
minimum- phase  transfer  functions.  A  practical  realization  of  the  algorithm  could 
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be  achieved  by  using  optoelectronic  picosecond  sampling  techniques  and  real-time 
digital  signed  processing.  Finally,  the  layer  impedance  profile,  using  measured  re- 


Figure  6.33  Characteristic  impedance  profile  for  multiple  discontinuities 

flection  data,  produced  excellent  impedance  estimates  for  both  the  initial  cascaded 
microstrip  step  discontinuity  and  the  narrow  0.3  cm  microstrip  line. 
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Summarizing,  we  have  presented  both  experimental  and  simulation  re¬ 
sults  for  three  microstrip  test  structures.  First,  pulse  propagation  simulations  were 
performed  using  both  the  equivalent  circuit  and  ARMA/AR  parametric  models  of 
an  IC  microstrip  interconnection.  Time-domain  measurements  were  compared  with 
the  simulation  results.  Second,  a  cascaded  microstrip  step  discontinuity  was  inves¬ 
tigated.  Using  the  equivalent  circuit  and  parametric  models  derived  in  Chapter  IV, 
reflection  and  transmission  response  simulations  were  compared  with  experimental 
results.  Third,  an  equivalent  circuit  model  was  developed  for  a  lossy  multi-section 
microstrip  structure  consisting  of  several  cascaded  step  discontinuities.  The  circuit 
model  was  validated  by  comparing  PSPICE  and  PUFF  simulations  with  microwave 
measurements.  Finally,  the  parameter  estimation  algorithms,  presented  in  Chapter 
V,  were  applied  to  simulated  and  measured  response  data  in  order  evaluate  their 
performance. 
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VII.  SUMMARY  AND  CONCLUSION 


With  the  development  of  automated  optoelectronic  measurement  systems  suit¬ 
able  for  the  manufacturing  environment  comes  the  problem  of  how  to  characterize 
microwave  integrated  circuit  (MIC)  devices  and  interconnections  in  MMIC’s.  This 
thesis  has  presented  new  computer-aided  design  models  and  parameter  estimation 
techniques  to  model  these  MIC  interconnections  using  picosecond  time-domain  mea¬ 
surements. 

A.  SUMMARY  OF  SIGNIFICANT  RESULTS 

Equivalent  lumped-distributed  circuit  models  were  derived  for  lossy  dispersive 
microstrip  transmission  line  structures.  The  circuit  models  consisted  of  resistive  and 
reactive  lumped-elements,  and  distributed  lossless  transmission-line  sections.  The 
models  were  compatible  with  general  computer-aided  design  techniques  and  simu¬ 
lation  programs,  such  as  PSPICE.  In  addition,  the  circuit  models  for  the  cascaded 
microstrip  step  discontinuity  accurately  characterized  the  reflected  and  transmitted 
wave  propagation  in  the  structure.  These  models  mark  a  significant  improvement 
over  empirical  modeling  for  which  accurate  generalized  closed-form  solutions  are 
not  available  for  CAD  programs. 

The  picosecond  sampling  performance  of  current  optoelectronic  measurement 
systems  motivates  a  digital  signal  processing  framework  for  the  modeling  and  anal¬ 
ysis  of  high-speed  transient  signals.  Parametric  models  of  lossy  dispersive  microstrip 
transmission-lines  provide  an  innovative  modeling  approach  not  previously  pre¬ 
sented  by  microwave  engineers.  The  parametric  model  approach  is  emphasized 
for  several  reasons.  First,  the  fitting  of  an  autoregressive  model  to  an  observable 
time  series  is  a  linear  process  which  can  be  handled  using  well  tried  and  highly 
efficient  computing  algorithms.  Secondly,  an  autoregressive  moving-average  model 
provides  the  smallest  number  of  parameters  for  an  optimal  and  most  parsimonious 
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representation  of  the  propagation  process.  Finally,  the  identification  problem  in 
parametric  modeling,  in  other  words  choosing  an  AR  or  ARMA  representation  for 
a  microstrip  line,  was  shown  to  be  intimately  related  to  the  critical  physical  prop¬ 
erties.  We  have  demonstrated  how  the  effective  microstrip  permittivity  of  a  MIC 
interconnection  can  be  evaluated  from  ARMA  or  AR  model  parameters.  Further¬ 
more,  parametric  model  simulations  accurately  characterized  the  picosecond  pulse 
propagation  on  lossy  dispersive  microstrip  structures.  A  first-order  AR  model  of 
the  cascaded  microstrip  step  discontinuity  formed  the  basis  for  the  development  of 
the  multi-section  microstrip  structure  model. 

Several  parameter  estimation  techniques  were  presented.  Both  deterministic 
and  stochastic  algorithms  were  developed  to  estimate  the  ARMA  and  AR  model 
parameters.  The  weighted  least  squares  (WLS)  algorithm  estimated  the  ARMA  dig¬ 
ital  filter  coefficients  or  model  parameters  from  a  prescribed  unit  impulse  response. 
The  algorithm  rapidly  converged  when  applied  to  several  finite-length  deterministic 
examples. 

A  discrete  transformation  of  the  network  function  approximation  method,  in¬ 
troduced  by  Elmore  [Ref.  46],  was  presented.  This  method  approximated  either  a 
first-order  “all  pole”  or  a  second-order  “pole- zero”  network  function  from  measured 
transient  pulse  delay  and  rise  times.  A  first-order  AR  model  and  a  second-order 
ARMA  model  were  derived  using  the  impulse  invariant  transformation  of  the  re¬ 
spective  network  functions. 

A  layer  reflection  coefficient  estimation  algorithm  was  derived.  This  algorithm 
estimated  the  reflection  coefficients  from  the  sampled  reflection  response  of  a  lossy 
multi-section  microstrip  structure.  The  algorithm  assumed  a  constant  attenuation 
factor  and  group  velocity  throughout  the  structure.  The  algorithm  detected  major 
discontinuities  within  the  structure. 

Finally,  the  layer-probing  algorithm  was  developed  to  characterize  multiple  dis¬ 
continuities  within  a  multi-section  microstrip  transmission-line,  consisting  of  several 


118 


step  discontinuities.  A  minimum-phase  parametric  model  described  each  discontinu¬ 
ity.  Inverse  filtering  and  deconvolution  techniques  were  employed  by  the  algorithm 
to  iteratively  probe  successive  sections  of  the  structure. 

B.  CONCLUSIONS 

The  research  conducted  for  this  thesis  has  shown  that  the  parametric  modeling 
approach  provides  an  alternative  framework  to  characterize  and  analyze  dispersive 
pulse  propagation  on  microwave  integrated  circuit  interconnections.  The  model¬ 
ing  of  dispersive  propagation  on  silicon  and  gallium  arsenide  microstrip  structures 
is  important  because  it  points  to  a  limitation  in  IC  speeds  even  when  transistor 
switching  speed  continue  to  improve.  The  measurements  conducted  and  reported 
in  this  work  have  verified  that  the  proposed  equivalent  circuit  models  are  compati¬ 
ble  with  available  computer-aided  design  (CAD)  programs.  The  circuit  models  were 
developed  using  a  combination  of  frequency-dependent  expressions  available  in  the 
literature  and  original  analysis  by  the  author. 

The  development  of  optoelectronic  measurement  systems  will  provide  the  ca¬ 
pability  to  directly  sample  high-speed  transients  on  the  semiconductor  substrates 
without  introducing  significant  parasitics.  With  the  availability  of  time-domain 
measurements,  we  developed  new  signal  processing  algorithms  to  estimate  both  the 
ARMA  and  AR  model  parameters. 

The  layer-probing  algorithm  presented  a  technique  to  iteratively  identify  indi¬ 
vidual  step  discontinuities  in  a  multi-section  transmission-line.  Simulation  results 
showed  how  destructive  interference  problems  caused  by  dispersion  and  multiple 
reflections  degraded  the  algorithm  for  some  structures  investigated.  This  problem 
was  overcome  by  measuring  and  modeling  the  transmitted  impulse  response. 
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C.  FUTURE  DIRECTIONS  FOR  RESEARCH 

The  objectives  of  the  research  were  successfully  accomplished.  Several  sugges¬ 
tions  for  future  research  are  now  presented. 

Researchers  are  currently  studying  optical  techniques  to  measure  nonlinear  and 
multiport  characteristics  of  MMIC  components  and  circuits.  Based  on  these  experi¬ 
ments,  nonlinear  parametric  models  should  also  be  investigated.  The  high  frequency 
sampling  capability  of  photoconductor  devices  can  be  integrated  with  pulsed  laser 
technology  to  develop  a  picosecond  time-domain  measurement  system  using  a  white 
noise  excitation.  Here,  several  stochastic-based  parameter  estimation  algorithms 
can  be  applied  to  the  modeling  problem. 

Equivalent  circuit  and  parametric  models  need  to  be  developed  for  coupled 
microstrip,  coplanar  waveguide  and  other  microstrip  discontinuities,  such  as  bends, 
crossings  and  the  asymmetrical  step  discontinuities. 

Improvements  to  the  layer-probing  algorithm  should  include:  1)  the  develop¬ 
ment  of  a  method  to  estimate  the  attenuation  factor  of  each  lossy  section.  2)  the 
implementation  of  a  minimum-phase  ARMA  parametric  model  for  each  section,  and 
3)  the  integration  of  cross-correlation  techniques. 

Another  important  aspect  of  this  research  is  its  extensions  to  other  engineering 
problems.  Future  efforts  envision  the  application  of  the  proposed  system  identifica¬ 
tion  algorithms  to  ultra- wideband  (UWB)  radar  target  identification.  Emphasis  is 
required  on  time- domain  electromagnetics  in  measuring  and  computing  the  impulse 
response  of  targets,  radar-absorbing  materials,  transmitting  and  receiving  antennas, 
and  various  other  media.  New  signal  processing  algorithms  must  be  developed  for 
impulsive  radar. 
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Appendix  A  Derivation  of  the  Schur  Recursion  Algorithm 

A  recursive  in  order  solution  for  the  AR  lattice  filter  reflection  coefficients  is 
presented.  This  derivation  of  the  so-called  Schur  algorithm  will  show  the  computa¬ 
tional  efficiency  of  this  algorithm  over  the  Levinson- Durbin  algorithm  [Ref.  12:  pp. 
6-22].  The  algorithm  derivation  will  assume  that  the  signal  data  is  stationary. 

1.  List  of  Symbols  and  Definitions. 

The  following  symbols  are  used  in  the  derivation: 
y(n)  The  autoregressive  (AR)  signal  sequence. 

b  The  prediction  error  filter  coefficient  vector. 

e(n)  The  forward  prediction  error. 

e(n)  The  backward  prediction  error. 

£{•}  Expectation  operator. 

[Rjj,]  The  autocorrelation  matrix  of  the  signal  y(n). 

B(z)  The  z-transform  of  the  forward  model  parameter. 

B(z)  The  z-transform  of  the  backward  model  parameter. 

The  (p  - f  l)-th  order  model  reflection  coefficient,  or  lattice  filter 
reflection  coefficient. 

E(z)  The  z-transform  of  the  forward  prediction  error. 

E(z)  The  z-transform  of  the  backward  prediction  error. 

The  (p)-th  order  normalized  forward  prediction  error. 
r*p*(?i)  The  (p)-th  order  normalized  backward  prediction  error. 
aSp'{n)  The  (p)-th  order  Schur  forward  recursion  parameter. 

0^(n)  The  (p)-th  order  Schur  backward  recursion  parameter. 
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2.  Algorithm  Derivation. 

The  present  output  data  y(n)  can  be  estimated  from  the  weighted  summation 
of  the  past  values  of  the  output  for  the  (p)-th  order  autoregressive  ( AR )  model.  This 
estimate  is  expressed  as  [Ref.  46:  p.  150] 

_  v 

yin )  =  YM  y{n  -  i)  (A. la) 

i=l 


or  in  matrix  form 

■'  rj~' 

y(n)  =  b  y(n) 


(A.  16) 


where 


bT  =  [6,  b2  ...  bp] 

yr(n)  =  [y(«  -  i)  y{”  -  2)  •••y(w-p)] 

The  prediction  error  is  defined  as  the  difference  in  the  signal  and  its  estimate  or 

e(n)  =  y(n)  -  y{n) 

=  y(n)  -  bT  y(n)  (  4.2) 

The  mean  squared  prediction  error,  expressed  as  a  function  of  the  prediction  error 


filter  weights  6,  is 

formulated  as 

[Ref.  56:  p.  20] 

S 

|e2(n)}  = 

bj  [ft,,  |  b  -  2  l/B.„ 

+  ^yy(O) 

where 

^yy(O) 

Ryyi  1)  . 

Ry  y(~P+  1  )' 

= 

Ryy(  1 ) 

«yy(0)  . 

Ryyip) 

.  RyyiP  —  1 

)  •  •  •  -Ryy(l ) 

Ryy(  0)  . 

and 

—yy  =  ^yy(^)  Ryyip)} 


The  solution  of  (A. 3)  for  the  optimum  filter  weights  is  given  by 

kept  ~  [-^-yy]  Qyy 
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which  is  commonly  call'''’  the  normal  equations  [Ref.  56:  p.  22],  The  normal 
equations  can  be  written  m  matrix  form  as 


Ryyi®)  Ryy(  1)  .  Ryyi  P"bl)  by  Ryy(^) 

Ryy{l)  ftyy(O)  .  Ryy(p)  ^2  Ryy(^) 


1)  Ryy{l)  Ryy(0)  J  |_  j 

The  autocorrelation  matrix  [Ryy]  is  Toeplitz  if  the  signals  are  stationary  which  is 
assumed  in  this  derivation.  The  Levinson  algorithm  provides  the  solution  of  the 
(p  +  1  )-th  order  model  parameters  by  using  the  previously  determined  p-th  order 
model  parameters.  [Ref.  46:  pp.  81-85].  The  (p  +  l)-th  order  model  parameters 
are  obtained  recursively  by 

B(p+X)(z )  =  B(p\z )  -  z~xk(p+l)B^p){z)  (A. 6a) 

B{p+X\z)  =  z'xB(p\z)  -  k(p+^B(p\z)  (A.66) 

where  B^p\z)  is  the  />  th  order  forward  model  parameter  transfer  function,  B^p\z) 
is  the  p-th  order  backward  model  parameter  transfer  function,  and  A-b’+1)  is  called 
the  reflection  coefficient.  The  z-transform  of  the  forward  and  backward  prediction 


error  is  obtained  by  multiplying  (A. 6)  by  Y(z). 

E[p+X){z)  =  E(p\z)  -  z~xk(p+x)>E{p\z)  (.4.7a) 

Eip+X){z)  =  z'xE(p\z)  -  k(p+x]E{p\z)  (A. 76) 

In  the  time  domain  the  prediction  error  equations  become 

e(p+1)(n)  =  e(p)(n)  _  k^p+x^e^-p\n  —  1)  (A. So) 

e(p+1)(n)  =  e(p)(n  -  1)  -  k^p+x)e{p\n)  (A.86) 


To  find  elp+x\n)  and  e^p+x\n)  from  e^(n)  and  eM(n),  respectively,  it  is  only 
necessary  to  determine  ^p+1).  Since  the  zero  order  prediction  of  a  signal  is  y{n)  =  0, 
the  zero  order  error  is 

e°(n)  =  e°(7?)  =  y(n)  (A. 9) 
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and  equations  (A. 8a)  and  (A. 8b)  can  be  realized  as  a  lattice  structure  as  shown  in 
Figure  A.l. 


Figure  A.l  Second  Order  AR  Lattice  Filter  Model 


The  reflection  coefficient  can  be  determined  by  minimizing  the  square  of  the 

expected  value  of  equation  (A. 8a)  and  (A.8b)  with  respect  to  the  forward  (kf)  and 
backward  (£*)  reflection  coefficients,  respectively.  The  forward  reflection  coefficient 
is 

s{eW{n)eW(n  -  1)} 


k/p+1)  = 

£  jeb^n  —  l)eW(n  —  1)  j 
and  the  backward  reflection  coefficient  is 

5{eW(n)eW(n  -  1)} 


(A.10) 


J bj(p+1)  = 


£|e^)(n)e(p)(n)| 


(All) 


The  forward  and  backward  reflection  coefficients  can  be  related  to  the  geometric 
mean  form  as  [Ref.  57:  pp.  434-437] 

5{eW(n)eW(n  -  1)} 


= 


^{[eWfn)]2  }£{[eW(n  —  l)]2} 


(.A.  12) 
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By  substituting  (A. 12)  into  (A. 10),  the  forward  reflection  coefficient  becomes 


(P+D  =  jAv^)\I£{  [e(P)(n  -  1)PT 


k  (p+i)  =  k(P+ 1) 

||  [c(p)(n  —  l)]2 


£{  [»)(»  -  l)f  ) 

e(p)(n)|| 


(,4.13a) 


(A.  13b) 


\eW(n)!\  =  j£{  [eW(n)]2} 


||e(p)(n  -  1)||  =  \JS{  [e(p)(n  -  1)]“  } 
Similarly,  the  backward  reflection  coefficient  (A.ll)  yields 

t (p+l)  = 

f{[eW(n)]2} 

tl(p+ 1)  = 

II  P*(»)f  II 


(A. 14a) 


(A. 146) 


Substituting  (A. 13b)  and  (A. 14b)  into  (A. 9a)  and  (A. 9b),  respectively  yields 


e(p+1)(n)  =  e^(n)  —  k^p+i^e^p\n  —  1)- 


||e<p)(n)|| 


(A. 15a) 


||  [eW(77  -  1)]' || 

e(p+1)(it)  =  e(p)(n  -  1)  -  fc(p+1)e(p)(n)  {A. 15b) 

ll[e(p)(«)]2|l 

The  normalized  (p  +  l)-th  order  forward  and  backward  prediction  errors  are  defined 


j(p+-)(n)  = 

1  ’  ||e(P+1)(n)| 


f(p+D(„  _  i)  =  1 


e(p+1)(n-  1) 


jje(f+1)(n  -  1)|| 

Factoring  ||e^(n)|j  from  the  right  side  of  (AA5a)  and  normalizing  yields 

i(p+'Hn] = ffSi|{8W(w)-  - »}  (■4-i7°i 


(A. 16a) 


(A. 166) 


Factoring  ||e(p)(rc  —  1)||  from  the  right  side  of  (A. 15b)  and  normalizing  yields 

||e(p)(n  — 


r^p+1^(n)  = 


||e(P+!)( 


n) 


,|f(p)(n  -  1)  -  fc(p+1)e(p)(n)|  (A. 176) 


We  will  next  show  that 


||e(p)(n)|| 

||e(P+1)(n)| 


1  - 


A'(p+i) 


and 


P(p)(n  ~  1)11 

||e(P+1)(n  -  1)|| 


(A.  18a) 


(A. 186) 


1  - 


fc(p+i) 


Proof:  Equation  (A. 18)  is  defined  as 


||£W(n)|| 


[f{  [eW(n)]2} 


||e(P+l)(n) 


(A. 19a) 


[£{[£(»+»(n-l)]2}]! 

Squaring  (A. 9a)  and  taking  the  expectation  yields 

£{  [e(p+1>(n)]2  }  =F{  [e<p)(n)]2}-21fc/p+1)5{e(p)(n)e^(n-l)} 

+  [fc(p+1)]2F{  [e(p)(n  -  l)]2  }  (A. 196) 

Substituting  (A. 10)  for  kf(p+}\  and  simplifing 

£:{  [e(p+1)(n)]2  }  =  £{  [e(pHn)]2  }  -  k  f{p+l)  8  [e(p)  [n)e{p)  (n  -  1)}  (A.19c) 

Returning  to  equation  (A. 19a),  we  substitute  (A. 19c)  to  yield 
||e(p)(n)||  1 


||e(P+1)(n) 


s{  [e(p)(n)]2  }  -  fc/p+1)£{eW(n)e(p)(n  -  1)} 
[e(P)(n)]2  } 


/fc/p+1)£{e(p)(n)e<p)(n  -  1)} 
£{  [e(p)(n)]2  } 


( A.19d) 
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Finally,  we  substitute  (A. 13a)  for  k/p+1)  into  (A.19d),  and  multiply  the  denomina¬ 
tor  term  si  le^(n)]  }  by  the  ratio 


yjs{  [e(p)(n)]2  [e<p)(n  -  l)]  '  } 


yj£{  eW(n)  2  e^(n  -  1)  2| 


simplifying  yields 


lMp)(n)il 

||e(P+1)(n)| 


||e(p+i)(n)| 


1 


1  - 


jt/(p+1)f|e(p)(n}e(;’)(n  -  1)J 

\jt{  [e(P)(n)]2  [ebd(n  -  l)]2  } 


1 

2 


1 


1  - 


k(P+ J) 


(A.19e) 


(A. 19/) 


Following  the  same  approach  using  equations  (A. 9b)  and  (A. 11),  it  can  be  shown 


that 


-  1)11 

||e(P+1)(n-  1)|| 


(A. 20) 


The  normalized  lattice  recursion  equations  are  from  (A.  17a)  and  (A.  17b) 


e(p+l)(n) 

1 

1 

_fc(p+i)‘ 

e^p\n) 

f(P+1)(n). 

y/l  -  [k(P+')}l 

— jfc(P+1) 

1 

_f(p)(n  —  1) 

or 


where 


e(p+l)(n) 

f(p+i)(n)  J 


e(p\n) 

[  f(pl(n  —  1 ) 


0[fc(p+1)]  = 


y/l  -  [k(p+  D]2 


1  -)t(p+1) 

_^(P+!)  1 


(A. 21a) 


(A. 216) 
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In  order  to  obtain  the  Schur  recursion  equations  we  will  introduce  two  new  param¬ 
eters,  a^(n)  and  p\n  —  1)  as 


«W(n)  =  f  |cW(Jb)  y(n)  j 


gSp\n) 


f  e(p)(n)  y(nh 
\  ||e(p)(n)||  / 

liiwif 


(A.22a) 

(A.22b) 

(A.22c) 


where  Ryy(n)  is  a  (p  x  1)  correlation  vector,  and  the  norm  ||y(n)|j  =  y/Ryy{Q). 


Similarly, 


fi(k)(p\n  —  1)  =  £  jr^(n  —  1)  y(n)  j 

_  5rr(p)(rc  - *)  y(n)i 

l  jjr(p)(n  —  1)||  i 


p(p),n  -  1)  =  _ - 

-  {  j  ||y(n-l)||* 


(A.23a) 

(A.236) 

(.4.23c) 


Using  the  normalized  lattice  recursion  equation  (A. 21)  and  the  previously  defined 
Schur  recursion  parameters  (a(p)(n))  and  (fftp\n  —  1)),  the  Schur  recursion  equation 
is 


a.(p+1)(n) 
/?(p+1)(rc  —  1) 


=  0(Jfc(p+1>) 


a(p)(n) 
£(p)(n-  1) 


(A.24) 


We  will  next  show  the  ((p+  1  ))-th  order  reflection  coefficient  in  terms  of  the 
ratio  of  Schur  recursion  parameters.  From  equation  (A. 22a)  and  (A. 23a)  we  form 


a(p>(n)  _  ^{e(p)(”)  y(n)} 

(3^p\n  —  1)  £|f(p)(n  —  1)  y(n)| 

£{e(p)(n)} 

£  jrO^n  —  I)  j 


(A. 25a) 

(A.256) 


128 


g{tW(n)}||8W(n-l)|| 

£{eW(„-l)}||eW(„)|| 


gSp\n) 

S' 

^p)(n)  e^(n  —  1)  j 

f 

-  1) 

V 

{mn)]2}^\ 

[  [e(p)(n 

-i)]2} 

(,4.25c) 


(A.25d) 


2  2 

Using  stationarity  and  |e^(n  —  1)J  J  =  |e^(n)j  J,  we  have  [Ref.  46:  p. 

168] 

aM(„)  £{eW(n)eW(n-l)} 

£(|,)(n-l)  £{[eW(n-l)]2j 


(A25/) 

(A.25g) 


3.  Normalized  Lattice  Filter  Realization. 

This  section  describes  the  normalized  lattice  filter  as  a  direct  method  of 
solving  given  the  autocorrelation  sequence  (Ryy)-  Assume  the  zero  order 

Schur  parameters  are 


Q(o>(i)  =  j?(i) 
0(o)(  0)  =  fl(0) 


Inserting  these  initial  values  in  (A. 25)  we  have 

*0)  _ 

_W) 

R(oy 

The  first  order  Schur  parameters  for  (n  =  1, 2)  are  solved  using  the  recursion  equa¬ 
tion  (A. 24)  as 

a(3)(l)  = _ 1 . .  [oW(l)  -  Ar(1)^°)(0)l 

yjl  -  I^1)]2  L  J 
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and 


<*W( 2)  =  -p=L=[gW(2)  -  Jfc(1)/3(°)(1)1 

y/i  ~  P>]2 1  J 

where  q^°^(2)  =  R( 2)  and  =  R(l).  The  /3^( n )  parameters  are 

0(1)(O)  =  — . =L==  f^(o)(0)  -  fc(1)a(°)(l)l 

vA-[*(1)]  J 


/3(1)(1)  =  -  =L=  f^°)(l)  -  fc(1)a(0)(2)l 

0  -  [fc(1)]2 

where  g^°^(2)  =  R{ 2)  and  /?(0)(1)  =  R{1).  Substituting  the  previous  results  into 
(A.25g)  yields 


q<1)(2) 

^>(0) 

R(2)-k^R{\) 


R(0)  -kWR(iy 


The  Schur  recursion  equation  can  be  realized  as  a  p-th  order  normalized  lattice 
structure  as  shown  in  Figure  A. 2. 


Figure  A. 2(a)  Zero  Order  (n  =  1)  Normalized  Lattice  Filter 


130 


{.  .  R(l),  R(0)) 


Figure  A. 2(b)  First  Order  (n  =  2)  Normalized  Lattice  Filter 
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APPENDIX  B  PSPICE  EQUIVALENT  CIRCUIT  LISTINGS 


IC  INTERCONNECTION  (PI-NETWORK  AND  MATCHED  IMPED  TERMINATION) 

• 

*  MICROSTRIP  TRANSMISSION  LINE  SIMULATION  (LCL1  PI  -  MODEL 

*  TEST  STRUCTURE  ( LENGTH  .  500E-06  m  ) 

*  SILICON  SUBSTRATE  ( ER  -  1 1 .7 ) 

* 

.OPT  ACCT  LIST  NODE  OPTS  NOPAGE  RELTOL-.OOI 

.WIDTH  IN-133  OUT-80 

.OPTIONS  LIMPTS  -  20000 

.OPTIONS  ITL5  -  50000 

.OPTIONS  TRTOL-5 

.OPTIONS  ITL4-5 

V  1  0  PWL(0  0  3.75p  .1 1  7.5p  .25  112Sp  1.45  15p  1.9  1 8.75p  2  22.5p  1 .9  26.25p  1 .85  37.5p  1 .6  52.5p  1 .2 
75p  .72  78.75 p  0) 

RS  1  2  74.23 

T1  2  0  3  0  ZO-74.23  TD-125ps 

R1  2  0  1.065K 

R?  2  3  10.4 

R3  3  0  1.065K 

Cl  3  0  0.03pF 

LI  3  4  0.332nH 

C2  4  0  0.03pF 

RL  4  0  74.23 

.TRAN  1.25pS  80ps  0  1.25pS 

.PROBE 

.END 

IC  INTERCONNECTION  (T-NETWORK  AND  MATCHED  IMPED  TERMINATION) 


.OPT  ACCT  LIST  NODE  OPTS  NOPAGE  RELTOL-.OOI 

.WIDTH  IN-133  OUT-80 

.OPTIONS  LIMPTS  -  20000 

.OPTIONS  ms  .  20000 

.OPTIONS  TRTOL-5 

.OPTIONS  m.4-5 

•  PIECEW1SE-UNEAR  APPROXIMATION  OF  NORMALIZED  PCE  SAMPLER  1 
V 1  0  PWL(0  0  3.75p  .11  7Sp  .25  1 125p  1.45  15p  1 .9  18.75p  2  22.5p  1.9  26.25p  1 .85  37.5p  1 .6  52.5p  1 .2 
75p  .72  78.75p  0) 

RS  1  2  74.23 
R1  2  01.065K 
R2  2  3  10.4 
R3  3  01.065K 

T1  2  0  3  0  20-7423  TD-223ps 
Cl  3  0  0.06pF 

T2  3  0  4  0  20-7423  TD-223ps 
RL  4  0  7423 

.TRAN  1 25  pS  1 0OpS  0  0.05ps 
.PROBE 

.PRINT  TRAN  V(4) 

.END 
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CASCADED-STEP  DISCONTINUITY  PSPICE  PI-MODEL  SIMULATION 
• 

*  DIELECTRIC:  (6-10  ER  -4.2) 

*  INPUT  SIGNAL:  (PIECEW1SE-LINEAR  APPROXIMATION  OF  THE  MEASURED 

*  GAUSSIAN  PULSE  OF  THE  HP851 0  NETWORK  ANALYZER) 

* 

OPT  ACCT 1 1ST  NODE  OPTS  NOPAGE  RELTOL-.OOI 

.WIDTH  IN-133  OUT-80 

.OPTIONS  UMPTS  -  20000 

.OPTIONS  ITL5  -  50000 

.OPTIONS  TRTOL-5 

.OPTIONS  ITL4-5 

V  1  0  PWL(0  0  100P  .06  300P  .17  400P  .37  800P  1.88  900P  2.02  INS  1.93  1.4NS  .41  1.5NS  .17  1.6NS 
.044 1.7NS  0) 

RS  1  24124 

T1  2  0  3  0  Z0-41.24  TD-750PS 
R1  3  03.118K 
R2  3  4 1.091 
R3  4  03.118K 

T2  4  0  5  0  ZO-41.24  TD-214PS 
Cl  5  0  7.13PF 
L4  5  6  0.362NH 
C2  6  0  7.13PF 

T3  6  0  7  0  20-41.24  TD-214PS 

R4  7  03.118K 

R57P1091 

R6  8  0  3.118K 

RL  8  041.24 

.TRAN  100PS  5NS  0  30PS 
.PROBE 

.PRINT  TRAN  V(2) 

.PRINT  TRAN  V(8) 

.END 
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CASCADED-STEP  DISCONTINUITY  PSPICE  SHUNT  CAPACITOR  MODEL  SIMULATION 

• 

*  DIELECTRIC:  ( EPOXY  G-10  ER  .  4.2 ) 

*  INPUT  SIGNAL-  (PIECEWISE-LINEAR  APPROXIMATION  OF  THE  MEASURED 

*  GAUSSIAN  HP 5510  NETWORK  ANALYZER  SOURCE) 

* 

.OPT  ACCT  LIST  NODE  OPTS  NOPAGE  RELTOL-.OOI 

■WIDTH  IN-133  OUT-BO 

.OPTIONS  UMPTS  -  20000 

.OPTIONS  ITL5  -  50000 

.OPTIONS  TRTOL-5 

.OPTIONS  m.4-5 

V 1  0PWL(0  0  100P.06  300P  .17400P  .37  800P  1.88  WOP  2.02  INS  1.93  1.4NS  .41  1.5NS  .17  1.6NS 
.044  1.7NS  0) 

RS  1  241.24 

T1  2  0  3  0  Z0-41.24  TD-750PS 
R1  3  0  2316K 
R2  3  4  1.409 
R3  4  0  2.31 6K 

T2  4  0  5  0  ZO-41.24  TD-214PS 
Cl  5  010.91  PF 

T4  5  0  6  0  ZO-7.424  TD-24.38PS 
T5  6  0  7  0  ZO-7.424  TD-24  .38PS 
C2  7  0 10.91  PF 

T3  7  0  8  0  ZO-41.24  TD-214PS 

R4802316K 

R5  8  9  1.409 

R6902316K 

RL  9  041.24 

•TRAN  100PS  5NS  0  30PS 
-PROBE 

.PRINT  TRAN  V(2) 

•PRINT  TRAN  V(9) 

.END 
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MULTI-SECTION  MICROSTRIP  FILTER  LAYER-PROBING  PSPICE  SIMULATION 
• 

*  DIELECTRIC-.EPOXY  (ER  -  4.3).  H-1.59mm.  DESIGN  FREQ  -  8.79GHz 

*  CHARACTERISTIC  IMPEDANCE  .  59.2  Ohms  (THEORETICAL) 

* 

.OPT  ACCT  LIST  NODE  OPTS  NOPAGE  RELTOL-.OOI 

.WIDTH  IN-133  OUT-80 

.OPTIONS  UMPTS  .  100000 

.OPTIONS  ITL5  -  100000 

.OPTIONS  TRTOL-5 

.OPTIONS  ITL4-5 

*  THE  FOLLOWING  PIECEWISE-LINEAR  APPROXIMATION  DESCRIBES  AN  INITIAL  100PS 

*  GAUSSIAN  PULSE  HAVING  A  3DB  BANDWIDTH  OF  8.79  GHZ 

*V3  1  0  PWL{0  0  10P  .03  20P  .33  SOP  .968  40P  1 .68  50P  2  60P  1.68  70P  .968  SOP  .33  90P  .03  1 0OP  0) 

* 

*  THE  FOLLOWING  PIECEWISE-LINEAR  APPROXIMATION  DESCRIBES  THE  SIMULATED 

*  INVERSE  FILTER  OUTPUT  USING  THE  1ST  DISCONTINUITY  AR  COEFFICIENT. 

*V3  1  0  PWL(0  0  10P  .03  20P  .302  30P  .668  40P  .8  50P  .47  60P  -.1426  70P  -.562  80P  -.552  90P  -.274 
100P  0) 

*  THE  FOLLOWING  PIECEWISE-LINEAR  APPROXIMATION  DESCRIBES  THE  SIMULATED 

*  2ND-ORDER  INVERSE  FILTER  OUTPUT  USING  U2  DISCONTINUITY  AR  COEFFICIENTS. 

*V3  1  0  PWL(0  0  1 0P  .03  20P  .276  30P  .398  40P  .206  SOP  -.238  60P  -.558  70P  -.436  80P  -.053  90P  .22 
100P  .212  110P  .024  120P  0) 

RSI  2 59.22 

TO  2  0  3  0  Z0-59.2  TD-1 P 
R1  3  0  831.7 
R2  3  4  8.48 
R3  4  0  831.7 

T2  4  0  5  0  ZO-59.22  TD-292  6PS 
Cl  5  0  3.57PF 

T4  5  0  6  0  ZO-59.22  TD-77.7PS 
R4  6  01.55K 
R5  6  7  4.52 
R6  7  01.55K 

T5  7  0  8  0  ZO-59.22  TD-77.7PS 
C2  8  0  1.5PF 

T6  8  0  9  0  ZO-59.22  TD-77.7PS 
R7  9  01.55K 
R8  9  10  4.24 
R910  01.55K 

T7 10  0  11  0  ZO-59.22  TD-77.7PS 
C3 11  0  3.57PF 

T8  1 1  0  12  0  ZO-59.22  TD-104.4PS 

R13  12  0  2.323K 

R14  12  133.02 

R15  13  0  2.323K 

RL  13  0  59.22 

.TRAN5PS23NS05PS 

.PROBE 

.PRINT  TRAN  V(2) 

.PRINT  TRAN  V(13) 

.END 
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